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I. INTRODUCTION 


Among the products of the collisions of 
the cosmic rays with the atmosphere which 
one might expect to find in detectable 
quantities in nature is tritium, radioactive 
hydrogen of mass 3, 12.5 half life, 18 years 
average life, which disintegrates to form the 
stable isotope of helium of mass 3. Tritium 
had been sought in natural waters prior to 
the discovery of its natural radioactivity. 
The methods used in this early search were 
inappropriate since they assumed that 
tritium would be a stable isotope and were 
not of sufficient sensitivity to detect natural 
tritium. It was fortunate, however, that the 
samples prepared for this original search of 
Lord Rutherford? and his colleagues at the 
Cavendish Laboratory were available for the 
search for radioactive tritium. It had been 
shown from earlier researches’ that the bom- 
bardment of nitrogen by fast neutrons of 
energies above 5 million volts produced 
tritium and the cross section for this reac- 
tion was such that the neutrons which had 
been observed by Korff* to be formed by 
the bombardment of the atmosphere by 
cosmic rays would produce tritium in suf- 
ficient yield to cause one to expect it to be 
observable.’ It was speculated at the time 
that the amount of tritium produced in this 
Manner could easily explain the extra- 
ordinary abundance of helium 3 in at- 


The Twenty-fourth Joseph Henry Lecture of 
the Philosophical Society of Washington, de- 
livered before the Society on March 25, 1955. 

*Lord Rutuerrorp, Nature 140: 303. 1937. 

*Cornog, R., and Lipsy, W. F., Phys. Rev. 
59: 1046. 1941. 

‘Rev. Mod. Phys. 11: 211. 1939. 

5’ Lippy, W. F., Phys. Rev. 69: 671. 1946. 
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mospheric helium. Helium 3 in atmospheric 
helium is about 1 part per million of the 
total helium, whereas in well or terrestrial 
helium it is only one-tenth as abundant as 
in atmospheric. Therefore, we understand 
that the helium 3 in atmospheric helium is of 
cosmic ray origin, probably. 

There is reason to believe that tritium 
produced by the cosmic rays would be pro- 
duced in the high levels of the atmosphere, 
probably above the tropopause, between 30 
and 50 thousand feet on the average. It, of 
course, could be produced at all levels to 
some extent though it would be most 
abundantly produced at the highest levels. 
We further expect that there would be op- 
portunity for most of the tritium atoms 
formed to burn to form water and conse- 
quently that the water of the atmosphere 
should be radioactive with tritium. An 
analogous research undertaken in Germany® 
at about the same time as our first work 
showed that the tritium might have some 
difficulty burning to form water by showing 
that atmospheric hydrogen is especially rich 
in tritium as compared to atmospheric water, 
it being about 10,000 times more concen- 
trated. However, it is still true that over 99 
percent of the tritium burns to form water 
since the atmospheric hydrogen is so ex- 
tremely rare relative to the moisture of the 
atmosphere. 

The Rutherford sample was obtained by 
Dr. A. V. Grosse and was electrolyzed and 
converted into deuterium gas, after which 
the deuterium gas was placed inside a geiger 
counter. It is necessary to count tritium in 


6 Fa.ttinas, V., and Harteck, P., Zeit. Natur- 
forsch. 6A: 438. 1950. 
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the internal volume of a geiger counter for 
the reason that its radiation is so extremely 
soft, 19 kilovolt energy upper limit to the 
beta spectrum. Because of this softness it 
will not penetrate with any reasonable ef- 
ficiency througha wall of any finite thickness 
nor emerge from a solid sample. The 
deuterium gas obtained by electrolysis and 
reaction of the heavy water produced with 
hot zine metal was mixed with a small pro- 
portion of ethylene gas (1 or 2 cm of mercury 
pressure) and a small amount of argon 
(about 3 em mercury pressure). This count- 
ing gas is a good geiger counter mixture and 
performs very well. The counter was shielded 
by standard anticoincjdence techniques and 
by a thick metal shield.’ The geiger counter 
had a natural background of about five 
counts per minute when so shielded though 
the volume was nearly eight-tenths of a 
litre. The electrolyzed heavy water sample 
produced for Lord Rutherford at an extreme 
enrichment of any tritium present of nearly 
ten million fold showed a very appreciable 
radioactivity. In fact, the radioactivity was 
so strong that it was several weeks before 
the failure of the counter to operate in a 
reliable fashion was understood by diluting 
the deuterium samples with ordinary hydro- 
gen gas and observing that the cause of the 
troubles was the excessive radioactivity of 
the deuterium gas from the Rutherford 
sample. The sample had been prepared by 
the Norsk Hydro heavy water concern by 
electrolyzing heavy water to further concen- 
trate any tritium which might have been 
present in the original water used to make 
the heavy water. The original water used to 
make the heavy water in the Rutherford 
sample came from the Lake Mésvan, which 
lies on the mountain plateau just east of 
Bergen in Norway. It was therefore quite 
certain that the water electrolyzed to form 
the Rutherford sample was fresh snow water 
and therefore had not had an opportunity to 
lose its tritium by radioactive decay. Dr. 
Grosse, Dr. Wolfgang, Dr. Johnston, and I 
all worried that the excessive radioactivity 
found in Rutherford’s sample might be in- 
advertent in that in 1936 when the sample 
prepared for Lord Rutherford’s search for 


7 ANpERSON, E. C., ARNoLD, J. R., and Lipsy, 
W. F., Rev. Sei. Instr. 22: 225. 1951. 
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tritium in nature had been studied in the 
Cavendish Laboratory heavy water was 
quite a rare commodity and it might have 
been used in the deuteron accelerator. 
Tritium is produced in prolific quantities by 
the reaction between deuterons and this 
tritium would have easily explained our 
observations. We therefore asked the Norsk 
Hydro Company to prepare a duplicate of 
the Rutherford sample. This was done, and 
it completely confirmed the assay found in 
the Rutherford sample. The abundance of 
tritium® using modern values® for the frae- 
tionation factor of the Norsk Hydro plant 
was 2.4 X 10-" atoms of tritium per atom 
of ordinary hydrogen, a value close to that 
expected? As a result of the discovery that 
rain and snow do, indeed, contain tritium in 
about the concentration expected, it was 
decided to pursue the assay of the natural 
waters of the earth for cosmic ray tritium 
and to study its possible practical useful 
applications. 


II. METHOD OF MEASUREMENT 


‘The water sample to be assayed for tritium 
content is first distilled in a standard 2-gallon 
gas-fired still, then is mixed with sodium 
hydroxide to produce a 3-percent (by weight) 
solution, and finally is electrolyzed in an 
iron-nickel anode plant'® to a volume of 
about 1 ce or less and the deuterium and 
tritium content of the final product deter- 
mined. [t has been observed that the ratio 
of the separation factors for tritium and 
deuterium can be taken as 2.0 within the 
experimental error of the data which is about 
5 percent." This results in the very simple 


expression, 
; ee of *. x Vo 
To No X V 


In this equation T and 7) are the final and 
original concentration of tritium respec- 





8 Grosse, A. V., Jounston, W. H., WoLrGana, 
R. L., and Lipsy, W. F. Science 113: 1. 1951. 
Material in this article was presented at the First 
Research Day of Temple University’s Research 
Institute, September 14, 1950. 

® Harteck, P., private communication. 

10 Brown, W. G., and Daacetrt, A. F. Journ. 
Chem. Phys. 3: 216. 1935. 

" KaurrMaNn, S., and Lipsy, W. F. Phys. Rev. 
93: 1337. 1954. 
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the ; = ; oe , j 
was tively in the water samples. N is the final yield considering that the intensity of the 
ave concentration of deuterium and No» is the cosmic rays themselves is somewhere be- 
tor initial concentration of deuterium which tween 0.5 and 1 primary particle per cm? 
a is 0.0156 percent for ocean water, Chicago ‘per second. In other words, between one- 
this | rain, 0.0149 percent, Chicago snow, 0.0141 third and one-seventh of all the primary 
our percent.’ If one takes the original volume of radiations hitting the earth succeed in pro- 
sk water to be 10-20 liters in the case of rains es a ae piginel cs high yields rn 
of | and rivers and 100 liters in the case of ocean Sheet a . D ed "Wolt es Poo 
. ‘ . rd C rR. L. Wolfge é 5 
and water and the final volume to be 0.3 to 0.5 “0YG Vurne, -T. : igs hpnendivate) 
i . . : sipinirai M. L. Kalkstein in progress at the present 
in ec. the resulting concentration of tritium, ; . ; 
: : time at the Brookhaven Cosmotron which 
of r, is adequate to allow several counts per . 9 0K hill; 
me ate alt ial bl cate cal th furnishes 2.05 billion volt protons and at the 
2) 2 OF ODS fable © , rate abdov . ° ay ae 
ant a 5 2 rai pincers a Se we Berkeley Bevatron which furnishes 5.7 bil- 
c a1 Fa € re : > . . . 
“i bac > equatonie: out our counts per minute. lion volt protons. The cross sections for the 
hat With this arrangement it has been possible production of tritium for a wide range of 
hat to complete a considerable number of as- targets encompassing the range of elements 
7 says for natural tritium.!4 rhe general result have been large. In fact, as an average result 
1 ee . : 
ai has been that the average tritium production one finds that about one-tenth of all the col- 
a rate by the cosmic rays is about 0.14 tritium lisions in this billion volt energy range 
wr: atoms per cm? per second of the earth’s sur- produces tritium. This is a large yield con- 
ful face, with a probable error of about 20 per- sidering the small binding energy which 
cent of this number. This is a very large tritium has and the essential instability of 
; ; ; al this particle. However, it does indicate that 
2 Craig, HARMON, private communication. ‘ , ; 
18 Figures quoted in earlier publication™ were the cosmic-ray production rate given by the 
incorrect. Numbers given here are the correct study of the tritium in the waters of the 
ye: members. earth is essentially correct. The assay results 
on 4 yonBuTTLaR, H., and Lipsy, W. F., Journ. ciceeneg ~ " eters , s 
sl Inorg. and Nuclear Chem. 1: 75. 1955. are given in Table 1. 
it) 
Ill. TABLE | 
an ARMS ake ae 
of Sample wer Tritium abundanc 
he — Description of sample (T idams por O° coon 
eT- 7: Chicago Rains and. Snows 
at) , eo oe hee P 
nd 5 May 11, 1951. c ollected from 1000 to 1200. Storm lasted from evening} 33 2 
h of 10th to afternoon of 11th. 3.81 in. rain. 
e 14 October 14, 1952. Collected throughout storm. Storm lasted from) 20.4 0.7 


ut 1700 to 2400. 0.70 in. rain. 
dle 16 November 17, 1952. Collected throughout storm. Storm lasted from! 37 
0030 to 0100; 0600 to 0630. 0.31 in. rain. 
17 November 18, 1952. Collected throughout storm. Storm lasted from) 66.0 
2100 of 17th to 1200 of 18th. 0.70 in. rain. 
18 November 22-24-25, 1952. Collected throughout storms. Rained| 
afternoon of 22nd, evening of 24th, all day 25th. 1.13 in. rain. | 
19 December 2, 1952. Fell during night. 0.29 in. snow. | -13.2 
January 6, 1953. Fell during afternoon. 0.05 in. snow. | ; 
22 January 23, 1953. Collected from 1515 to 1530. Storm lasted from| 9.0 
1200 to 1900. 0.32 in. rain. 
30 February 11, 1953. Collected from 0930 to 1245. Light rain all inn 5.0 
1. ing. 0.03 in. rain. | 
; February 16, 1953. Collected 1030. Fell from 0300 to 1400. 0.19 in. 
nt snow. 
34 | February 20, 1953. Collected 1615. Storm lasted from 1530 to 1645. 3.3 
0.94 in. rain. 
38 March 3, 1953. Collected throughout storm. Rain, sleet, and snow} 9.4 
fell from 1300 to 2000. 0.38 in. rain. 
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no. 
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132 


135 
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TABLE 1 (C ‘ontinued) 





Description of sample 


Tritium abundance 
(T atoms per 10'* H atoms) 





March 7, 1953. Collected 1200. Fell from 0800 to 1500. 0.18 in. anow.| 

March 12, 1953. Collected 1030. Storm lasted from 0900 to 1700. 0.98 
in. rain. 

March 14, 1953. Collected throughout storm. Rained from 1000 to! 
1600; 2000 to sees. 1.05 in. rain. 

March 18, 1953. Collected throughout storm. Light rain from 1300 to! 
1700. 0.13 in. rain. - 

March 21-22, 1953. Collected throughout storms. Rained night of| 
21st, night of 22nd. 0.06 in. rain. 

March 31, 1953. Collected from 1030 to 1100. Rained lightly from 
1000 to 1200. 0.07 in. rain. Wtd. Ave. (December to March): | 

April 3, 1953. Collected from 1000 to 1200. 0.07 in. rain. 

April 15, 1953. Collected from heavy rain 1130. 0.60 in. rain. | 

April 30, 1953. Rained heavily from 1830 to 1915 with thunder -" 
lightning. 0.87 in. rain. 





| May 22, 1953. Collected at 0945. Heavy thunderstorm. 0.82 in. rain. | 
| June 5, 1953. Collected from 1530 to 1610. Thunderstorm. 0.43 in, rain.) 


June 25, 1953. Collected from intermittent rain between 2130 and| 
2400. 0.02 in. rain. 
July 1-2, 1953. Collected between 2300 on July 1 and 0300 on July 2.| 
Trace rain on July 1, 0.06 in. on July 2. | 
July 5, 1953. Collected from July 3 to July 6. Storm occurred about! 
2200 on July 5. 0.39 in. rain. | 


| July 17-20, 1953. Collected from 0930 July 17 to 1600 July 20. Rain-| 








| December 2, 1953. Rain. Collected at 1500. 


| January 21, 1954. Snow. Collected at 1700. 


storms on morning of July 17, afternoon of July 17, night of July! 
18, and morning of July 20. Rainfalls were ;,0.57 in. July 17; trace 
July 18; 0.10 in. July 20. 

August 3, 1953. Collected 1930 to 1945. Trace rain. | 

August 4, 1953. Collected 1445 to 1455. 0.06 in. rain. 

September 4, 1953. Collected 0400 to 1100. 0.50 in. rain. This was the! 
first rain after 3 weeks drought. 

September 18, 1953. Rained very hard—thunderstorm. Collected at 
2145. 0.63 in. rain. 

October 18, 1953. Trace rain. 

October 26, 1953. First drops collected. Total rainfall was 0.31 in., 
but this sample consisted of the first drops. Collected at 1530. 

October 26, 1953. Same storm as sample No. 102 except water col- 
lected after three hours of continuous rain. Collected at 1800. 0.31 
in. rain. 

November 20, 1953, Collected 9000 to 1100 from the first part of the 
rainstorm. 0.56 in. rain. 

November 20, 1953. Collected 1600 at end of rainstorm. 

November 27, 1953. Collected in the morning after an overnight 
snow. 


December 2, 1953. Rain. Collected at 1700. 
December 3, 1953. Rain. Collected from 0925 to 0930. 
December 12, 1953. Rain. 

January 20, 1954. Rain. Collected from 0845 to 0855. 





January 26, 1954. Rain. Collected from 1000 to 1400. 

February 5, 1954. Rain and Snow. Collected after 1700. 

February 15, 1954. Heavy rain with thunder and lightning. Collected 
from 1900 to 2200. 

February 16, 1954. Snowy rain. Collected from 1500 to 1600. 

March 2, 1954. Snow fell 1500 to 2400. Shoveled from the ground on 
March 3, 1954 at 1430. 

February 20, 1954. Rain. Collected from 1230 to 1315. 

March 19, 1954. Rain. Collected from 1430 to 1540. 
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TABLE 1 (Continued) 


Tritium abundance 


Sample Description of sample (T atoms per 10'* H atoms) 


ms) no 


167. | March 24-25, 1954. Rain. | 283 
168 March 29, 1954. Fresh snow. 196 
169 March 29, 1954. Collected March 31, 1954, after most of snow had| 145 
melted. 
170 April 7, 1954. Heavy shower. Collected 1530. | 248 
171 April 15, 1954. Short, heavy storm. Rain. Collected 1130-1145. | 360 
174 April 19, 1954. Collected at 1745. It had been raining since about! 425 
noon. | 
175 April 20, 1954. Collected at 0900 during rain on April 21, 1954. | 260 
181 May 26-27, 1954. Collected during rainfall at night. | 450 
182. | May 27, 1954. Coilected at 1105-1110. Different rainstorm from at 416 
ple No. 181. 
183 June 1, 1954. Collected 0930-0935. Rained previous night and all day.| 390 
Average tritium content from December 1, 1952 to December 1, 1953) 8.2 
weighted according to the total rainfall for stcrm samples for 27| 
samples amounting to 12 inches out of the total rainfall of 26) 
inches. 
Unweighted average for the samples between December 1, 1952, and| 7.8 
December 1, 1953. 
Calculated average for rainfall over the past 18 years (Lake Michi- 7.4 
gan). 


or we we 


t HHH HHH HEE 





B. ‘Rains and Snows from Other Areas 





42 enetvetite, hein: snow. January 23, 1953. 1.17 in. snow. 5.5 

4 Wilson Springs, Arkansas, rain. Collected on February 10, 1953 in a| 2.2 
heavy thunderstorm. 

49 Honolulu, Hawaii, rain. Collected on morning of March 26, 1953 “ 0.61 
Honolulu, 5329 Keikilani Cirele by Mrs. Terry Yoshida. 

95 Hickham Field, Oahu, Hawaii, rain. Collected at Hickham Field, | 0.59 
July 17, 1953. 

98 Manila, Philippine Islands, rain. Collected 5 June, 1953 “ 
the Bureau of Quarantine by Dr. R. Abriol. 

110 Tantalus, Koolau Range, Oahu, Hawaii, rain. Collected 13 Novem-| 0.73 
ber 1953. The rain was collected by residents of Tantalus. | 

120 Oahu, Hawaii, rain. Collected by pooling the catches of several rain| 0.67 
gages located at: Lower Laukaka, Nuuahu Reservoir No. 4, Wil.) 
Rise No. 4, Manoa Valley, Palolo Valley, Nuuahu Pali, Kalihi, 
U.S8.G.S. Station, Kaliha, Tunnel No. 2, Kaliha Reservoir Site, 
Tantalus Peak, Bureau of Water Supply, City and County of| 
Honolulu. 

121 Manoa Valley, Hawaii, rain. December 1953. Single rain. 

125 Wellington, New Zealand, rain. Collected October 1, 1953 on the 
roof of the American Embassy at Wellington, New Zealand. The} 
roof had recently been painted and consisted of galvanized iron.| 
The meteorology: a vigorous depression was located at the Tas-| 
man Sea west of Cook’s Strait. It moved ESE and came over| 
Wellington in the early hours of October 2, 1953. Rainfall from the} 
forward side of the depression began when the cloud base was at| 
3 kilometres. After it began to rain the cloud base moved down to} 
1 to 2 kilometres and the freezing level was at 2 kiometres./ 
The clouds were alto stratus and nimbo stratus. 

126 Harwell, England, rain. Two gallons collected after a windless os 22.5 + 0.6 
on October 24-25, 1953. Artificial tritium from Harwell may have) 
been present in this sample. Sample sent by J. M. Fletcher. 

127 Harwell, England, rain. Collected October 26 to 28, 1953. Artificial) 25.5 + 0.4 
tritium from Harwell may have been present in this sample. | 

128 Santa Barbara, California, rain. Rain fell on January 23 and Jan-| 4.0 + 0.1 
uary 24, 1954. Collected by Dr. V. L. Vanderhoof on January 24,| 
1954 from Mission Creek, Santa Barbara. Total rainfall was 3) 
inches. 
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TABLE | (Continued) — 
i Sa 
ie Descspsten of seaphe en aT... , 
164 Kobe, Japan, rains. November, 1953 (4 or 5 storms). Sent by Prof. 6.5 + 0.4 
Seishi Kikuchi, Osaka University. 
165 Puerto Rico, rain. March 26, 1954. Drained off a tin roof at Ramez 65.6 + 1.2 
Air Force Base. Total rainfall 0.42 inches. 1 
179 Valparaiso, Chile, rain. April 4 or 5, 1954. After 3 months’ drought. 4.3 + 0.2 
Rained very heavily for about 36 hours. : 
Water from Lake Mésvan used to supply the Norsk-Hydro. The 2.4 1 
Lake Mésvan water probably was melted snow from the winter of 
1946 and 1947. The water was taken from the lake at the end of L 
January, 1948. (8) 1 
zz L: 
C. Mississippi River 
= li 
23 Mississippi River, St. Louis, January 31, 1953. 5.6 + 0.6 
24 Mississippi River, Rock Island, January 29, 1953... | 2.5 + 0.3 — 
26 Mississippi River, St. Louis, February 4, 1953. 4.5 + 0.6 
27 Mississippi River, Rock Island, February 6, 1953. 3.7 + 0.4 — 
28 Mississippi River, Memphis, February 4, 1953. 6.0 + 1.0 
29 Mississippi River, New Orleans, February 8, 1953. 4.7 + 0.3 
31 Mississippi River, St. Louis, February 10, 1953. 6.0 + 0.9 
36 Mississippi River, St. Louis, February 20, 1953. 6.4 + 0.5 
37 Mississippi River, Rock Island, February 24, 1953. 4.4 + 0.2 
44 Mississippi River water from Rock Island. Collected March 16, 1953. 3.2 + 0.2 1 
47 Mississippi River water from St. Louis. Coliected 0900 on March 17, 5.4 + 2.4 9 
1953. 
57 Mississippi River water from Rock Island. Collected 1315 on April 5.3 + 0.3 9 
17, 1953. 10 
58 Mississippi River water from St. Louis. Collected 1300 on April 17, 6.0 + 0.4 
1953. 
80 Mississippi River water from Rock Island. Collected on June 30, 7.2 + 0.7 
1953. 14 
88 Mississippi River water from St. Louis. Collected 1320 on. July 22, 7.3 + 0.4 
1953. ‘ 
Average for Rock Island 4.7 14 
Average for St. Louis 6.0 
D. Other Rivers 14 
10 Columbia University distilled water. Collected August 5, 1944. 2.4 + 0.1 
(4.5 + 0.2 as of 1953) —_ 
15 Sangamon River, Decatur, August 6, 1952. 1.15 + 0.08 
52 Arkansas River, Conway, Arkansas, March 20, 1953. 3.12 + 0.10 aS 
104 Riber Elbe. Collected at Hamburg, Germany, on August 31, 1953. 2.57 + 0.12 5 
105 River Weser. Collected at Bremen on September 1, 1953. 1.76 + 0.10 
108 River Rhone. Collected near Lyon, France, on September 10, 1953. 2.64 + 0.16 9 
111 River Main near Wiirzburg, Germany. September 13, 1953. Collected 1.76 + 0.19 
by Dr. H. V. Buttlar, who collected all the European samples. 
112 River Loire. Collected at Digoin, France, on September 9, 1953. 2.11 + 0.14 15 
114 English stream near river Cam. Water taken about one mile from the 1.25 + 0.10 
source (a spring). Collected near Cambridge, England. 15 
122 River Donau. Collected near Ulm, Germany, on September 12, 1953. 2.13 + 0.38 
123 River Mosel. Collected near Metz on September 7, 1953. 2.15 + 0.12 15 
124 River Seine. Collected near Nogent, France, on September 8, 1953. 1.80 + 0.3 16 
131 River Fulda. Collected near Kassel on September 24, 1953. 2.35 + 0.1 16 
132 River Rhine. Collected between Geisenheim and Riidesheim on 3.0 + 0.3 
September 7, 1953. 
133 River Marne at Joinville, collected September 8, 1953. 2.1 + 0.2 
143 Shasta Dam, California. January 30, 1954. 2.7 + 0.1 
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no 
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40 


143 


148 
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53 


93 


159 
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161 


TABLE 1 (Continued) 


Description of sample 


EK! Rito de los Frijoles, Jemez Mountains, New Mexico. Collected 


February 7, 1954, at Bandelier National Monument. Sample prob- 
ably represents average over the winter snowfall in the Jamez 
Mountains. Sent by IX. C. Anderson. 

Rio Grande, northwest of Santa Fe, February 7, 1954. The spring 
thaw was beginning to raise the river level. but it had not yet 
reached peak volume by any means. Sent by E. C. Anderson. 

Winsor Creek, taken February 22, 1954, at Cowles, New Mexico, just 
above junction with Pecos River. Sent by E. C. Anderson. 

Rio Guajataca, Puerto Rico, at Lares, collected March 2, 1954. 

Rio Arecibo, Puerto Rico, at Utuado, collected March 2, 1954. 

River Tomokoa, Florida, on Route 92 near Daytona Beach, collected 
March 19, 1954. 

Alafia River, Florida, on Route 60 about 20 miles east of Tampa, 
March 22, 1954. 


EK. Lake Michigan and Other Lakes 

Distilled water from the Jones Chemical Laboratory at the Uni- 
versity of Chicago. The intake for the water supply for the 
southern part of the city is two miles off the shore in southern 
Lake Michigan. Collected on May 12, 1952. 

Jones Laboratory, Tap water, July 7, 1952. 

Oak Park tap, hot water heater, ca 12 years old. 

Jones Chemical Laboratory distilled water. Collected on February 
6, 1953. 

Jones Laboratory, Tap water, February 13-16, 1953. 

Jones Chemical Laboratory tap water. Collected on October 26, 1953. 

Average for Lake Michigan 

Calculated average rain from the depth and area of the lake and the 
mean assay for tritium in the lake: 

Lake behind Shasta Dam, California, collected by Dr. L. R. Libby 
probably January 30 or 31, 1954, certainly not later than February 
15, 1954. 

Roundout Reservoir, south shore, collected February 6, 1954, by 
Seth Harris, Lamont Observatory, 500 yards in front of gate off 
New York 55. Ice cap over whole reservoir. 

Shandakan Tunnel effluent from Schoharie Reservoir at Allaben, 


New York, right off route N.Y. 28. Collected February 6, 1954,| 


by Seth Harris. 


F. Hot Springs 


Water from the main reservoir at Hot Springs, Arkansas. Collected 
on March 18, 1953. 

Water from the Lardarello, Italy, hot springs near Pisa. These hot 
springs furnish about % of the electric power inthe whole of Italy. 
The water actually was steam from the volcanic fumaroles. 

Bowers Hot Springs, Bowers Mansion, Nevada, collected on March 
11, 1954. 

South Steamboat Well, Steamboat Springs, Nevada. Collected 
March 11, 1954. 

Spring No. 24, Steamboat Springs, Nevada, collected March 11, 1954. 

Spring No. 50, Steamboat Springs, Nevada, collected March 11, 1954. 

Wilbur Springs, Colusa County, California, collected probably on 
March 23, 1954, by Donald E. White, T = 127.2°F. This water is 
exceptionally high in alkali chlorides, bicarbonates, boron, sul- 
phide, iodine, and other compounds. 
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TABLE 1 (Continued) 





Description of sample 





163 
176 


177 
185 


186 


41 
99 
89 


and 
119 


62 
63 


64 


69 


Devil’s Kitchen and Range Spring, The Geysers, collected March 


23, 1954, by D. E. White. T = 202°F. No surface discharge, but 
probably slight sub-surface. Vigorous boiling, evaporation from 
pool about 15 feet in area. Acid reaction. 

Magnesia Spring, The Geysers, California, collected March 23, 1954, 
by D. E. White. T = 122°F. Discharge estimated 2 gpm. 

Hot Springs, Arkansas, main reservoir. Collected at Government 
free Bath House, April 29, 1954. 

Wilson Springs, Arkansas, drill hole No. 1, collected April 29, 1954. 

Morgan Spring, Lassen Park. Growler Spring. This is near neutral 
alkali chloride-boron-bicarbonate-sulphate spring opposite in 
type from sample 186. T = 205.4°F; discharge about 8 gpm. Water 
believed to be in part voleanic but dominantly of surface origin, 
circulating deeply and mixing. Collected June 8, 1954. 


Sulfur Works, Lassen Park. Big boiling spring. This is an acid spring} 


typical of most of the Lassen hot springs. T = 190°F, no dis- 
charge. These springs appear to result from partial condensation 
of steam that is probably meteoric in origin. Some direct con- 
tribution of recently fallen snow. Collected June 5, 1954. 
G. Well Water 

Drill core water from the Wilson uranium prospecting area in Wilson 
Springs, Arkansas. Taken from the drill core on March 6, 1953. 

Water from the 900 feet deep well at the McDonald Observatory at 
Fort Davis, Texas. Sample given by Dr. Gerard Kuiper of Yerkes 
was pumped in June of 1953. f 

Champaign-Urbana city well water. Water taken in July 1953 from 
well No. 51, located west of Champaign. This well is 296 feet deep 


and is producing at a rate of 2,225 gallons per minute. The first) 


hundred feet of formation is Wisconsin glacial drift and the next 
200 feet are Illinoisan and Kansan drift. The very bottom of the 
well at 300 feet is at the top of the Pennsylvatiia system. No- 
where in the entire log of the well is any clay formation which 
would prevent the downseepage of rain. It was therefore thought 
that perhaps the well could contain some rainwater. It is gen- 
erally thought that the well water which supplies the area is 
ancient and possibly of melted glacier origin. Sample furnished 
by T. E. Larson, Head, Chemistry Sub-division of State Water 
Survey. 

H. Vintage Wines 


Widmer’s New York Riesling wine, Vintage 1952, Naples, New York 
Widmer’s New York Riesling wine, Vintage 1940, Naples, New York 
Widmer’s New York Riesling wine, Vintage 1946, Naples, New York 
Hermitage Rhéne wine, Vintage 1929, Tain, Dréme, France. 
Hermitage Rhéne wine, Vintage 1942. Tain, Dréme, France 
Hermitage Rhéne wine, Vintage 1917, Tain, Dréme, France. 
Hermitage Rhéne wine, Vintage 1951, Tain, Dréme, France. 


Chateau Laujac Bordeaux wine, Vintage 1928, France. 


Tritium abundance 
(T atoms per 10! H atoms 


10.8 + 0.4 


2.06 + 0.15 
0.57 + 0.88 


0.17 + 0.11 


—0.48 + 0.37 
(age greater than 50 
years ) 

0.13 + 0.05 
(age greater than 5 
years) 


5.3 + 0.3 

(5.6 + 0.3) 
3.2 + 0.2 

(6.6 + 0.4) 
3.63 + 0.16 

(5.4 + 0.3) 
1.13 + 0.38 

(4.3 + 1.4) 
2.15 + 0.21 

(3.92 + 0.4) 
2.15 + 0.28 

(3.0 + 0.3) 
3.4 + 0.4 

(3.8 + 0.5) 
1.16 + 0.16 

(4.6 + 0.7) 
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TABLE 1 (Concluded) 
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Semple Description of Sample a en 
70 Chateau Laujac Bordeaux wine, Vintage 1934, France. 1.16 + 0.30 
(3.3 + 0.9) 
71 Chateau Laujac Bordeaux wine, Vintage 1939, France. 2.6 + 0.4 
(5.6 + 0.9) 
72 | Chateau Laujac Bordeaux wine, Vintage 1945, France. 2.70 + 0.18 
(4.2 + 0.3) 
83 Sherry, Vintage 1942, Jerez de la Frontera, Spain. 1.93 + 0.27 
(3.55 + 0.5) 
84 Sherry, Vintage 1947, Jerez de la Frontera, Spain. 1.99 + 0.55 
(2.67 + 0.75) 
85 Sherry, Vintage 1951, Jerez de la Frontera, Spain. 2.73 + 1.0 
(3.0 + 1.0) 
K. Sea Water Samples 
75 Sea water collected from the beach at Santa Monica, California, on 0.54 + 0.02 
June 8, 1953. 
115 Sea water collected from the surface of the Gulf Stream on Septem-; 0.19 + 0.05 
ber 11, 1953, at 38° 5’ N, 69° 30’ W. The water temperature was 
82°F. This was collected by the ship, Vena, by Mr. Bruce Heezen, 
Chief Scientist from the Geophysics Section of the Geology De- 
partment, Columbia University, the Lamont Laboratory. This 
is Lamont Sample No. T-32. 
116 Atlantic sea water taken in Sargasso Sea at 34° 00’ N, 62° 50’ W. 0.29 + 0.02 
Temperature of the water was 81.5°F. Collected following a light 
drizzle. Lamont Sample No. T-28. 
117 Atlantic sea water taken from 300 miles from New York on the con- 0.55 + 0.04 
tinental slope at 30° 20’ N and 71° 30’ W. The water temperature 
was 75°F. Lamont Sample No. T-34. 
154 Atlantic sea water from the Sargasso Sea at 34° 00’ N. 52° 35’ W. 16 + 0.1 
Temperature of the water was 81.5°F. Collected September 7,| 
1953, immediately following a light drizzle of 1 hour. Lamont 
Sample No. T-25. 
L. Miscellaneous 
9 Cistern, Decatur, collected August 6, 1952. Covered loosely about 29) 5.9 + 0.5 
years ago. 
12 Cistern, Sullivan, Illinois, collected August 6, 1952. Covered with) 2.9 + 0.2 
tight iron lid about 14 years ago. (6.3 + 0.4) 
20 Fire extinguisher, Skokie, Illinois. Filled June 5, 1936. 12.5 + 0.5 
(32 + 1.2) 
7 0.54 + 0.02 


5 Pacific Ocean, Santa Monica, California. June 8, 1953. 





IV. DISCUSSION 


Recently it has been shown that the 
tritium produced by the cosmic rays is pro- 
duced only in part by the neutron reaction 
mentioned earlier,> a large fraction of the 
tritium being produced by the direct inter- 
action of the primary cosmic rays with the 
air.'° It would seem reasonable that whatever 
the mechanism be by which the cosmic rays 
produce tritium by bombarding the air, the 


Phys. Rev. 91: 922. 1953. 


1’ FiREMAN, E. L. 





production rate should vary with latitude in 
a manner not too dissimilar from that in 
which the secondary neutrons vary.'* Simp- 
son has published data for the variation in 
neutron intensity with latitude at 30,000 
feet altitude. Using these data, we can calcu- 
late the expected ratio of the worldwide 
average production rate of tritium, Q (T 
atoms per cm? per second) to the local pro- 
duction rate, Q, for various localities. For 


16 Simpson, J. A., Jr., Phys. Rev. 88: 1175. 
1951; 84: 335. 1951. 

















example, at Lake Mésvan, Norway, the 
local Q should be multiplied by 0.58 to 
obtain Q. At the approximate center of the 
Mississippi Valley the factor would be 0.64, 
and at Chicago it would again be 0.58. In 
this way the data given in Table 1 can be 
treated and corrected for the expected varia- 
tion in latitude, though this variation has 
not yet been completely proved for natural 
tritium. 

The most striking new development in the 
nature of the results of the assay of natural 
tritium is included in section K of Table 1, 
the data on surface sea-water samples. 
Whereas we had originally supposed that the 
waves would certainly Mix the sea to such an 
extent that the tritium would be undetect- 
able in surface sea water, we find that indeed 
this is not so and that surface sea water does 
contain tritium. Four results given in Table 
1, the beach water at Santa Monica (sample 
75), the continental shelf Atlantic water 
(sample 117), the Gulf Stream sample 
(sample 115), and the Sargasso Sea (sampie 
116) all strongly indicate that there is 
definitely measurable tritium in the surface 
waters. In order to understand this we should 
know that uniform mixing to a depth of 113 
meters would give an average T-value just 
equal to Q. This result is obtained if one 
assumes that all tritium finds its way into 
the sea eventually before decaying into 
helium 3—in other words, that the storage 
in ground water and in the atmosphere is 
negligible relative to the runoff and direct 
rain into the oceans. Since, as we will show 
later, there is little doubt that Q is between 
0.1 and 0.2 T atoms per cm? per second, one 
observes that the results in Section K of 
Table 1 would correspond to a mixing depth 
of about 100 metres. This comment assumes 
that the two coastal water samples are high 
because of rains from the neighboring conti- 
nents, which are richer, and takes the 
samples from the Gulf Stream and the Sar- 
gasso Sea as being more typical of the open 
ocean. It seems likely that © is probably 
less than 0.25, the average of the Gulf 
Stream and the Sargasso Sea samples, so 
we may expect that future measurements 
will show that the sea mixes to a depth of 
about 50 metres in the lifetime of tritium. 
This result is in agreement with the notion 
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that some stratification exists in the sea 
below the surface at a depth of 50 to 100 
meters. Measurements of the isotopic oxy- 
gen in air dissolved in Pacific Ocean water 
as a function of depth showed a point of 
inflection in the isotope ratio at a depth of 
about 100 meters,” and the Baltic is known 
to behave as almost completely stratified 
with the bottom of the surface layer lying at 
about forty meters.'* 

Though further measurements will be re- 
quired to establish the tritium content of 
the surface of the ocean, it is clear that this 
tritium content is low, and water vapor 
evaporated from the sea is_ essentially 
tritium-free. It probably would be well to 
calculate Q assuming sea-water vapor on 
the average to have about 0.25 T atoms per 
10 H atoms. Using this we can set up an 
equation which expresses the balance be- 
tween the production and decay processes: 


7, (T, — 0.25) = 4.7Q, (2) 


where r, is the average annual precipitation 
over the seas in meters, T, is the average 
tritium content of sea rain in T atoms per 
10° H atoms, and Q is the cosmic-ray tritium 
production rate in T atoms per cm? per 
second. The factor 4.7 is a numerical con- 
stant arising from the units. Equation (2) 
simply assumes that no appreciable storage 
of tritium in the atmosphere for times com- 
parable to the lifetime of tritium—18 years, 
average life—occurs, and that essentially all 
the tritium forms water. The latter assump- 
tion is well justified by the measurements of 
the tritium contents of atmospheric hydro- 
gen that have been made, which indicate 
well less than 1 percent of the tritium is 
found in this form. The value of f, used is 
0.77 meter per year.’ Using the data in 
Table 1 and the latitudinal variation factor 
described above, we see that T, should lie 
between 0.9 and 1.0. This gives Q values of 
0.11 and 0.12, respectively. 

This calculation involves the implicit 
assumption that the tritium made over the 

17 RaKEsTRAW, N. M., Rupp, E. P., and Mr. 
Doe, Journ. Amer. Chem. Soc. 73: 2976. 1951. 

18 GusTAFSON, T. and KULLENBERG, B. Svenska 
Hydrografisk-Biologiska Kommissionens Skrifter 
Ny Serie: Hydrografi XIII. 


® RaNKAMA, K. and Sanama, Th. G., Geo- 
chemistry, Chicago, 1950. 
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oceans is precipitated into the oceans. We 
know, however, that on the western coasts 
of the continents considerable tritium from 
the skies over the seas precipitates on land. 
Also, our measurements show that rains on 
the east coast of continents are much higher 
in tritium. Thus the sea rain just east of the 
eastern continental coasts must be much 
richer in tritium than T, as though land- 
born tritium were precipitating into the 
sea. In view of the prevailing westerly direc- 
tions of the winds in the regions measured, 
these results suggest that Q may be some- 
what higher than indicated by the ocean 
data. 

We can make similar calculation for the 
land areas by using the Mississippi Valley 
average of 6 multiplied by the latitudinal 
factor 0.6 as a worldwide average for land 
rain. 

It seems that the continental rains must 
be richer than sea and coastal rains, for two 
reasons. First, the length of time that the 
moisture is exposed to contamination by 
cosmic-ray tritium after leaving the ocean 
is greater, and secondly, there is less moisture 
in the air over the continents than over the 
oceans and coastal regions. This raises the 
specific activity, that is, the number of 
tritium atoms per hydrogen atom. 

We thus obtain 

Bi(6 X 0.6 — T,) = 4.7Q (3) 
where Pp; is the worldwide average runoff 
from the land, and the average composition 
of ocean rain is subtracted to correct for the 
tritium content of the sea-water vapor com- 
ing in over the western coasts and the Gulf 
of Mexico. The value of p,; used is 0.28."* The 
runoff must be used in equation (3) because 
re-evaporation allows a given. tritium atom 
to be precipitated more than once over land. 
Only transport into the sea where wave ac- 
tion dilutes the tritium accomplishes the 
fixation on the surface, and even in this case 
as we have seen previously the fixation is not 
complete. The value of Q so derived is 0:16. 
Therefore we conclude Q is likely to be close 
to 0.14, probably to within 20 percent. 

Using this result, we calculate the ex- 
pected tritium contents for average ocean 
rain and surface-ocean waters on the as- 
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These results are given in Table 2. It is 
pleasant to observe the agreement with the 
data given in Table 1. 


TABLE 2.—CaLcuLATED LocaL Tritium Propuc- 
TION Rates, Q, AND ExpEcTED OCEAN 
RAIN AND SurRFACE WATER TRITIUM CONTENTS 


$1 3| 2 
=| £) & 
ae 
Localit > | ag | ts 
Locality 2) ae| £2 
, E | 28) 3: 
=| sz] 35 
=) =| ge] && 
sic ie |& 
CNN on. 5h eicsi sl cats .. -/0.58'0.24 |2.0 10.55 
Norway and the North At- : 
lantic Ocean...... 0.55/0.26 |2.2 (0.58 


Mississippi Valley....... .. .|0.60/0.23 |1.9 |0.53 

Honolulu and the Hawaiia | 
OS oe eee 

New Orleans, San Francisco, 
Southern Italy, Southern 
Spain.... 

Puerto Rico...... 

Marshall Islands. 


_ 


-90/0.073}0.65,0. 17 


.00/0.14 |1.2 |0.32 
.25/0.17 |0.96\0.26 
.40/0.058)0.48/0.13 


a) 





It is obvious that the data in Table 1 for 
samples collected later than March 15, 1954, 
are not germane to our principal point—the 
cosmic ray tritium. These high numbers 
were the consequence of the thermonuclear 
tests in the Pacific in the spring of 1954. 
Most of our considerations in this paper will 
ignore them, except in connection with the 
samples from the hot springs. In this case, 
the principal question, whether the springs 
run rain water or not, can be answered nearly 
as well with the contaminated rain as with 
the normal rain. In fact, the appearance of 
unusually radioactive water in the effluent 
water in several springs must indicate turn- 
over times of a few days or weeks, at most, 
since the test series began on about March 1. 
The one rain from Chile collected on April 
4, 1954, showed no contamination, while all 
Northern Hemisphere rain and river waters 
tested showed it for March 15 onward. Ap- 
parently the mixing of water vapor across the 
equator is much slower than is the east- 
west mixing. 

The large fluctuations in the tritium con- 
tents of successive rains in the Chicago area, 
excluding the large rise in March 1954 due 
to the test activities in the Pacific, require 
special consideration in a meteorological 
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sense. Our basic equations (2) and (3) for the 
relations between the tritium content, T, and 
the cosmic-ray production rate, Q, are based 
on the assumption that tritium cannot stay 
in the air for years even if it is formed well 
above the tropopause—the stagnant layer at 
30,000 to 50,000 feet between the strato- 
sphere and the troposphere in which the 
normal weather phenomena occur. Half of 
the tritium probably is produced on the 
average above the tropopause. We assume 
that it cannot stay up there for times com- 
parable to 18 years—the average life of 
tritium. This assumption does not mean, 
however, that it may not mix horizontally 
in an east and west andin a north and south 
direction rather well due to stratospheric 
winds before penetrating through the tropo- 
pause and being precipitated. Whatever evi- 
dence for latitudinal variation of Q is to be 
found in Tables 1 and 2 should be taken as 
indicating that the stratospheric storage 
time is less than the north and south horizon- 
tal mixing time. 

Assuming perfect vertical mixing, we can 
speak of the atmospheric moisture content, 
w, in meters of water per cm’, and the 
storage time, 7 in years, for the moisture in 
the air, counting total time elasped in the 
air after the moisture leaves the ocean and 
until it returns eventually to the sea as rain 
or snow or runoff river water. These quanti- 
ties then will be related to the tritium con- 
tent of rain, T,, and the production rate Q, 
for the open sea by 


%- % = — (4) 


where Ty is the tritium content of ocean- 
water vapour. For the continental areas this 
simple equation does not hold, because re- 
evaporation of precipitated moisture without 
dilution is an important phenomenon. 

The values of w are not well known, but 
it seems likely from available studies?®: *! 


20 Jacoss, W. C., The energy exchange between 
sea and atmosphere Bull. Scripps Inst. Oceanogr. 
6(2): 22-122. 1951. 

21 Benton, G.S., Estoque, M. A., and Domin- 
11z, J., Johns Hopkins Univ. Department of Civil 
Engineering Scientific Report no. 1, Contract 
AF 19(122)-365 of the Geophysics Research 
Division of Air Force Cambridge Research Center. 
1953. 
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that the range of values will be 0.01 to 0.3, 
with corresponding average 7 values of 
0.013 and 0.039 or 4.7 days and 14 days, 
respectively. At an average wind velocity 
of ten miles per hour these numbers corre- 
spond to transport distances of 1,100 and 
3,300 miles respectively. 

It is clear that in general land rain 
should be richer in tritium than sea rain 
because the atmospheric residence time 
should be longer and because the mositure 
content, w, should in general be lower. 
However, it remains to apply the equa- 
tion to air masses of known w and with 
known trajectories so that the variation of 
Q with latitude can be introduced. It is 
to be hoped that careful coordination of 
meteorological data with tritium assays 
will yield results of importance to meteor- 
ology in general. 

An attempt has been made to understand 
the variations of the tritium activity in 
single rainstorms. Forty samples of indi- 
vidual storms were collected on the campus 
of the University of Chicago over a period of 
more than a year. These samples correspond 
to 16.23 inches of rain, or 46 percent of the 
total rainfall during this period. The T 
average of these rains, weighted by the 
amount of rainfall, gives 8.8 units, which 
compares reasonably well with the Lake 
Michigan assay of 7.4. (The Lake Michigan 
assay of 1.64 is used to calculate the com- 
position of the rain supplying the lake.) 

Approximate values for w were computed 
from radio sounding data for the period of 
December 2, 1952, to November 20, 1953: 
monthly means from the Weather Bureau 
Technical Paper no. 10 had to be employed 
for the period from November 20, 1953, to 
March 2, 1954. Values for w(T — To) are 
given in Table 3. 

While the measured values for T show a 
minimum in winter, it can be seen that this 
effect is over-compensated by multiplica- 
tion with the amount of precipitable water 
in the air, thus yielding the total tritium ac- 
cumulated in the air mass. One then ob- 
tains a maximum in summer. This is prob- 
ably due_to the following causes: (1) In 
summer the wind velocities are lower, there- 
fore the air mass is exposed longer to the 
tritium production. (2) It is known that the 
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continent is a water donor during the sum- 
mer-time. This means that the air mass 
picks up tritiated water from the continent. 
In winter the large water sources for the 
atmosphere is the ocean, which is almost 
dead in tritium. 


TABLE 3.—RAINFALL DATA AT THE 
UNIVERSITY OF CHICAGO CAMPUS 


| 
| Monthly 









































Date T | eee pe i | Joliet, ill w (T —To) 

| s/em* | "for w in | 

| or w in 

g/cm? 
12-2-62| 13.2 | 12.95] 1.2 | 1.0 | 15.5 
16-53 | 7.1 | 6.85| 0.5 | 0.8 | 3.4 
1-23 | 9.0 | 8.75| 1.8 | 15.7 
2-11 5.0 | 4.75| 2.0 | | 9.5 
216 | 10.8 | 10.55| 0.7 | 0.85 | 7.4 
22 | 3.3 | 3.05| 2.1 | 6.4 
$3 | 9.4 | 9.15| 1.0 9.2 
37 6 | 9.6 | 9.35| 0.8 | 7.5 
$12 | 4.75| 4.5 | 2.2 | 0.98 | 9.9 
$18 | 7.9 | 7.65| 2.1 | 16.1 

| i 
3-21/22| 3.1 | 2.85| — } = 
331 | 9.5 | 9.25| 2.5 | | 23.1 
= 6| 62 |. aes) 18°} | 15.9 
415 |10.0 | 9.75| 1.5 | 1.48 | 14.6 
42% | 9.0 | 8.75| 4.2 | | 36.7 
430 | 6.6 | 6.35| 2.6 | 16.5 
5-22 6.7 | 6.45| 4.8 | 2.15 | 30.9 
65 | 6.18| 5.9 | 3.25] 3.2 | 19.2 
625 | 7.9 | 7.65) 4.1 | 31.4 
720 | 4.3 | 4.05) 4.5 | 18.2 
75 | 7.6 | 7.35| 3.9 | 3.44 | 28.6 
7-17/20| 7.1 | 6.85) — -- 
8-3 3.6 | 3.35] 4.9 | 3.3 | 16.4 
8-4 17.0 | 16.75 | 4.1 68.8 
4 | 7.3 | 7.05| 4.3 30.3 
—| ——— 2:5 | 
918 | 12.0 | 11.75 | 3.3 38.8 
10-18 | 9.8 | 9.55] 1.9 | 1.8 | 18.2 
10-26 | 19.9 | 19.65 | 2.6 51.1 
11-20 | 7.3 | 7.05| 2.8 20.4 
11-20 | 13.5 | 13.25 | 1.2 | 15.9 
11-27 | 34.5 | 34.25 | 41.2 
12-2/3 | 7.3 | 7.05 1.0 7.1 
1-20-54 | 10.5 | 10.25 | 8.2 
1-21 | 18.7 | 18.45 | 0.8 | 14.7 
1-26 «=| 19.4 | 19.15 | 15.3 
a ae a 
2-5 | 23.0 | 22.75 | 19.3 
215 | 9.5 | 9.25) | 0.85 | 7.9 
216 =| 15.6 | 15.35 | 13.0 
220 | 42 | 3.95) | 3.4 
3-2 | 20.8 | 20.55 | 0.98 | 20.1 
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The argument that the vertical mixing 
may not be complete would yield a high 
value in summer also, since thunderstorms in 
summer are known to reach higher up than 
snowstorms in winter. But our measure- 
ments do not support this thesis, since 
thunderstorms do not average higher than 
rains, while snow in general seems to be 
quite active. 

The data given in Table 1, section H, on 
the tritium contents of vintage wines, to- 
gether with the comparison of the Lake 
Michigan average as shown in section E 
with the general average for the rainfall in 
the Chicago area, show that without doubt 
the tritium content of rain up to March 15, 
1954, was essentially the same as it had been 
for the last eighteen years. It is clear, also, 
of course, that if the calculation can be made 
for Lake Michigan in which the hydrology is 
known, and if the checks are satisfactory, 
that in the case of an unknown lake the 
measurement of the tritium and a compari- 
son with the rainfall for the general region 
should give the hydrology, or at least some 
data on it. For example, the average depth 
might be calculated, or the average storage 
time for water in the lake. 

The application of natural tritium to the 
determination of the storage or holdup 
time for ground waters may be of consider- 
able importance. Some applications have 
been made and are given in sections F and G 
of Table 1. In section F data are given for 
natural hot springs. It seems likely that all 
of the hot springs utilize surface water except 
the Lardarello (sample 93) near Pisa, which 
actually was steam from volcanic fumaroles, 
the Wilbur Springs at Colusa, Calif. (sam- 
ple 161), and the Morgan Spring at Lassen 
Park, Calif. (sample 185). 

The wells tested and reported in section 
G gave old water in agreement with ex- 
pectations. The lower limit of the age of 
these waters is 50 years. Thus we cannot 
distinguish water older than 50 years from 
very old or juvenile water. 

It is clear from an examination of the 
data in section H of Table 1 that tritium 
allows one to determine the age of wine. It 
is also clear that the water in the wines has 
essentially the same tritium content as the 








rain of the general area. For example, the 
Widmer Winery near Rochester, N. Y., 
shows a general tritium assay of 5.8 + 0.3, 
in good agreement with what one might 
expect in view of the results for the Missis- 
sippi Valley as a whole. The French wines 
in the Rhéne Valley near Bordeaux run 
around 314 or 4 of our units; Spanish wines 
about 3, which within the experimental 
error agrees with the general rainfall average 
for the area which is shown by the river- 
water assay. Therefore it is relatively cer- 
tain that agricultural products can be dated 
and also that the presence of rainwater in 
the agricultural products can be determined. 
It would be interesting, of course, to test 
for pumped irrigation water which might 
possibly be ancient in character, just how 
effectively the water was incorporated into 
the crop. 

To balance the cosmic-ray production of 
tritium, there must be some mechanism for 
the product of the radioactive disintegra- 
tion of tritium, helium 3, to escape from the 
earth. Being helium, no chemical compounds 
can be formed, and being nonradioactive, 
no nuclear transformations. Therefore, every 
second for each square centimeter of the 
sarth’s surface, 0.14 helium 3 atoms must 
escape on the average if the cosmic rays 
have been constant in intensity. A direct 
analysis of the air shows that there are 1.5 X 
10" helium 3 atoms per cm? of the earth’s 
surface. Therefore, we calculate that on the 
average a helium 3 atom must stay on the 
earth 34 million years. This escape time 
may well be about half this, since it is 
probable that the cosmic rays can also pro- 
duce helium 3 directly as well as through 
the tritium, but it seems unlikely that this 
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would do more than double the total rate 
of production of helium 3. Therefore we 
rather definitely can say that if the evs- 
mic rays have remained constant for the 
past 50 or 100 million years, the amount of 
helium 3 on earth is such as to indicate an 
escape time of between 20 and 40 million 
years. This is a rather reasonable magnitude 
considering the probable escape time for 
helium 4, as judged from other evidence. 
One finds in the acceptability of this result 
some evidence that the cosmic rays indeed 
have remained constant for this great pe- 
riod of time. 

It is by no means obvious that on a time- 
scale of 18 years the cosmic rays should re- 
main completely constant, if any appre- 
ciable portion of therm are connected with 
the solar sunspot cycle of 11 years’ duration. 
The data given in this report so far, however, 
indicate that there is no real evidence that 
the present rate of cosmic-ray generation of 
tritium is any different from what it has 
been over the past 18 years. 
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MATHEMATICS.—Spectroscopic eigenvalue analysis. C. Lanzcos,! Dublin In- 
stitute for Advanced Studies, Dublin, Eire. (Communicated by R. K. Cook.) 


In honor of Lyman J. Briggs on his eightieth birthday 
(Received May 20, 1955) 


A method is herein described which trans- 
forms the search for the real eigenvalues 
and eigenvectors of a matrix into the search 
for hidden periodicities of a function com- 
posed of periodic components. The applica- 
tion of the Fourier transform displays the 
entire eigenvalue spectrum of the matrix 
in the form of sharp maxima of a Fourier 
spectrum, in analogy to the operation of a 
spectroscope. The method is well suited for 
the big electronic calculators and has the 
advantage of simple programming, com- 
bined with high precision. The same method 
yields an iterative solution of large-scale 
linear algebraic systems. 

1. The search for hidden periodicities. The 
spectroscope is a physical instrument con- 
structed for the purpose of resolving an 
arbitrary superposition of periodic functions 
into its components. Let f(t) be of the fol- 
lowing form: 


f(t) = a cos yt + +--+ + Gy, COS Upt 
(1.1) 
+ b, sin vt + --- + 5b, sin v,¢. 


Then the spectroscope transforms this func- 
tion of ¢ into a new function F(v) of the 
frequency v, in the following manner. At 
certain definite values y = v; sharp lines ap- 
pear, the “spectral lines,’ whose intensity is 
proportional to +/a? + b?. In view of the 
finite resolving power of the instrument, 
however, the lines are not sharp but have 
a finite width. This means that the light 
intensity is not concentrated at the discrete 
frequencies vy = »; but falls off continuously 
according to the law 


sin (v — v,) rN 
(v — vi) #N 


The larger the N, the greater is the resolving 
power of the spectroscope. 

The following method of analyzing the 
eigenvalues and eigenvectors of a symmetric 


! Former staff mathematician, National Bureau 
of Standards. 


matrix (or more generally of any matrix 
whose eigenvalues ‘are all real) is called 
“spectroscopic” because it imitates mathe- 
matically the operation of a spectroscope. 
The great accuracy of spectroscopic ob- 
servations is caused by the high resolving 
power (large N) of spectroscopic instru- 
ments. In our operations the value of N will 
remain within limits which are modest in 
comparison to the values encountered in 
physical spectroscopy. And yet the high 
precision of spectroscopic measurements will 
be realizable (order of magnitude 107°), 
because of the great accuracy with which 
the basic function f(t) is available. 

Let us assume that f(t) is given at the 
equidistant values of ¢. 


(1.2) t = 0, 7, 27, --- , Nr. 
We introduce the notation 
(1.3) v, = 6, 


and restrict all 6; to the range [0, x]. More- 
over, ior our present purposes the sine 
analysis will not be needed. Hence we assume 
that the following sequence of ordinates is 
at our disposal: 


Ye = a cos kO,; + a cos kb, 
(1.4) + --- +a, cos k@, . 
d=@6L2 +> Ww) 


The aim of the spectroscopic analysis (also 
called “‘search for hidden periodicities’) 
may be characterized as follows: Given the 
ordinates y» = Yo, Yi, °°* , Yn, find the 
hidden frequencies »% (or the 6 which are 
proportional to them), and the associated 
amplitudes a; . 

2. Solution by Fourier analysis. The solu- 
tion of our problem can be given by the 
method of the Fourier transform (ef. [1], 
p. 99) suitably modified for the case of dis- 
crete data. The process of integration has 
thus to be changed to a process of swmma- 
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tion. We define the following function of the 
continuous variable é: 


P(e) = 4 yo + yr cos © § 


(2.1) 
+ Ye cos 5 2+:---+lyw cos 5 NE. 


This function is periodic with the period 2N. 
Moreover, it is an even function of €: 


F(—£) = F@). 


Hence it suffices to reduce the range of & 
to [0, N]. We will restrict ourselves to 
integer values of §: =k, which are of par- 
ticular interest. Hence we will put 


F(k) = u = 4 wo 
(2.2) Be 


‘a > Ya COS a v k +l yw cos N +k. 


These « can be plotted as isolated ampli- 
tudes at the equidistant abscissas k = 
0, 1, 2, --- , N. The sharp maxima of these 
amplitudes? will reveal the frequencies 
which are present in the original function 
f(t) and the strength with which they are 
represented. 

We shall first write down the Fourier 
transform F(¢) of the specific function 
yx = a; cos ké;. However, for the purposes 
of mathematical operations it will be con- 
venient to replace the frequencies @; by 
proportional quantities, to be denoted by 


Di: 
(2.3) 6, = P. N ° 


Let us introduce the following function of x: 


(24) (x) = 4 sin xx cot SH 


Then we obtain 


(2.5) F(t) = o(p — &) + o(p + §). 


2 For the sake of distinction we will refer to 
the yx (equidistant ordinates of the given function 
f(t)) as the “‘ordinates,’’ and the u (equidistant 
ordinates of the Fourier transform) as the ‘‘ampli- 
tudes.”’ 
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For all practical purposes the function (2.4) 
can be replaced by the simpler function 


(26) o2) - Tes 
ax 





An ideal situation exists when all the p; of 
a given problem are integers. Then 


us 
2.7 6; =m;=, 
(2.7) N 
where m; is an integer between 0 and N. 
The Fourier transform u, of the function 


(2.8) - Ya = 4; COS ™M; v a 


is now reduced to a single line at the point 
k = m;, with the amplitude xa . The dis- 


tribution of the u consists now of isolated 

lines. The position and magnitude of these 

lines determine directly the 6; and the a; . 
In the general case we can put 


(2.9) p=m+e., 


where m is an integer and ¢ a quantity be- 
tween 0 and 1. In this case we no longer ob- 
tain a single line for the Fourier transform 
of a cosine function, but a group of lines 
clustered around the maximum at k = m, 
with slowly decreasing amplitudes. We will 
write 


(2.10) k=m+u 


(u = integer) and obtain as the transform 
of the function cos 6a the amplitudes 


N sin ex (—1)* 
Um+u = DY —<—< «ee 


(2.11) 
w €—yp 
The sharp peak at » = 0 is still there but 
the neighboring amplitudes are no longer 
zero. They fall off rather slowly (if ¢ is near 
to 14), with alternating signs, and influence 
the position of other maxima, if these max- 
ima are not separated by large distances, 
ie., if N is not excessively large. 

We can greatly reduce the mutual inter- 
ference of neighboring peaks if instead of 
operating with the amplitudes themselves we 
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investigate each combination of three neigh- 
boring amplitudes. The second difference of 
the function 1/(e — 2) falls off with the 
third power of « — x and hence becomes 
quickly negligible. In view of the factor 
(—1)* in (2.11), the second difference has 
to be replaced by the second sum: 


(2.12) vu = Usa + 2m + Ung. 


The sign of the uw generally alternates be- 
tween + and —, but this regular + pattern 
is occasionally disturbed by a + +, or — — 
sequence. This irregularity reveals the near- 
ness of a peak, regardless of the magnitude 
of the u,. We underline each + + or — — 
pair and draw all further conclusions from 
them and their neighbors. The order m of 
the amplitude u,,, belonging to the first 
member of the pair, constitutes the integer 
part of p. The fractional part ¢ is now deter- 
mined as follows. We can put 
Mu 
(2.13) Unin = ( -1) 
et. 


Then 


“aa 


aa 28 
— e(1 — (2 —)* 


We now form the ratio 


= Vis 2-— e 
(2.15) q = api 0 
Um+1 l a € 
whence 
. = 
(2.16) ¢ = 7 r q. 
q 


The strength 8 of the peak is evaluated as 
follows: 


(2.17) B= (1 — ©)(Um + Um4i). 


Finally, going back to the original frequency 
and amplitude of the cosine component of 
the function (1.1), we get 
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2: 
(2.18) “* 92? 
and 
(2.19) 6; = (m+ 6 y 


This method of resolving a function, given 
in equidistant ordinates, into its periodic 
components is more accurate than the 
customary numerical schemes pursued in 
the “search for hidden periodicities’ (cf. 
[9], p. 349). It is shaped to the demands of 
a problem in which the basic ordinates y. 
are given with high accuracy; we need not 
rely on the assumption that the frequencies 
which are present in f(t) are necessarily 
widely separated. If N is sufficiently large 
(order of magnitude 1,000) and the func- 
tional ordinates given with an accuracy of 
10°, we can obtain an accuracy of 10° or 
higher, for the determination of the 6; . 

3. The Chebyshev polynemials. Let x be 
an ordinary algebraic quantity, limited to 
a real positive range between 0 and 1. We 
generate a succession of polynomials by the 
following recurrence relations (cf. [10], p. 
XI): 


bo = 1, bh = 1 — 22 
by = 2(1 — 2x)bi — bo 
(3.1) b; = 2(1 — 2x)be — dy 


by = 2(1 = 2x)by-1 — by-2 


The quantities thus generated have the fol- 
lowing significance: 


(3.2) b. = cos ké, 
where 
(3.3) sin’ 5 = 2. 


We now replace z by the matrix A. If 
we assume that the eigenvalues of A are all 
real, positive, and limited to the range 
{0, 1], we may construct the “operating 
matrix” 
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(3.4) C =21 —4A 


and generate a succession of vectors by the 
following iteration scheme; where bp is an 
arbitrary trial vector: 


b, = 14Cbo 
bo = Ch. — bo 


(3.5) b, = Ch, — by 


by - Chy-1 —_ by-2 ‘ 


Let us analyze bo in the reference system of 
the principal axes of -4. We will call these 
principal axes w, , We, *** , Wn; 


(3.6) be = Biwi + Bow. + +--+ + Brn. 
Then 

( b. = (Bi cos k0,)w, ae (Be cos ke) we 
3.7 


+ --- + (8, cos k0,)w,, 


where the @; are in the following relation to 
the eigenvalues \,; of the matrix A: 

(3.8) A; = sin’ - 
We will now pick out one definite component 
b,—e.g., the first component of each suc- 
cessive vector b,—and consider this one- 
dimensional sequence of quantities as the 
ye-sequence (1.4): 

(3.9) ye = bf”, bf”, B®, --- , o& 

We subject this sequence to a Fourier analy- 
sis, obtaining the amplitudes (2.2). The 
isolated maxima of this sequence, evaluated 
according to sec. 2, will determine all the 
6—and thus also all the A;—with high 
precision. 

4. Programming for the large electronic 
calculators. In contrast to the author’s 
earlier method of ‘minimized iterations” 
(cf. [4]; see also [2] and [8]), the present 
method is far less economical in the number 
of iterations employed. The total number of 
iterations N must exceed the order n of the 
matrix by a factor of 10 to 20, in order to 
guarantee the proper independence (and 
thus easy separability) of the peaks. On the 
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other hand, the present method recomme ids 
itself by the simplicity of the operations and 
the ease of its programming for the elec- 
tronic computer. Since the successive itera- 
tions follow each other quite automatically 
without any disturbance, the generation of 
the b, vectors is a quick process, even in the 
case of large matrices. Moreover, in the 
previous method the eigenvalues \,; did not 
appear explicitly but had to be obtained as 
the roots of an algebraic equation. In the 
present procedure, the eigenvalues come 
directly in evidence by the position of the 
maxima of a set of Fourier amplitudes. 
Furthermore, the usual difficulty of properly 
protecting the small eigenvalues against the 
encroachment on the part of the large 
eigenvalues is here completely circumvented 
by the method of projecting them on a 
semicircle with the radius 144. On the circle 
the question of magnitude loses its meaning, 
and all eigenvalues come into play on an 
equal footing. In this manner, a single 
Fourier analysis is able to display the entire 
spectrum of a matrix. 





OP i oe 
A = 4 — 16 cos 0 


The procedure consists of two phases: 

1. The generation of the vectors b, . 

2. The Fourier analysis of these vectors. 
The latter analysis need not be carried out 
with the entire vector but only with a 
single component of each of these vectors. 
In actual practice, it is advisable to isolate 
two sets of components. We can use the 
second set for checking purposes and also 
for the testing of multiple eigenvalues, as 
we will see later, in sec. 9. 

The symmetry of the matrix A is not 
required. It suffices to know that all the 
eigenvalues of A are real.’ However, the 


_ * The case of complex eigenvalue will be treated 
in a subsequent paper. 





Oc 


sy! 
am . 
are 
fine 
ne\ 
thi 
the 
casi 
bey 
nee 
cap 
dec 
nec 
ber: 
con 
V1 
fror 
the 
the 
5. 
orig 
exce 
eige 
estil 
Ger: 
abs¢ 
choc 
num 
the | 


(5.1, 


we ( 
mat) 

If. 
all t! 


(5.2) 


Gers 
the 

disac 
sque: 
semi 
used. 
a pr 
evalt 
eigen 
preci 
value 
not 1 
are | 


. 10 


- ids 
and 
lee. 
era- 
‘ally 
n of 
the 
the 
not 
1d as 
the 
ome 
the 
des. 
erly 
the 
arge 
ited 
na 
ircle 
‘ing, 
) an 
ngle 
ittre 


out 
h a 
ors. 
late 
the 
also 


not 
the 
the 


ated 





OcTOBER 1955 Cc. 


symmetry of A has an important advantage. 
If A is symmetric, then the eigenvectors w; 
are orthogonal to each other. Since we de- 
fined the vectors b, by a process which 
never multiplies the eigenvectors by any- 
thing larger than 1, we know in advance that 
the absolute value of the vectors b, (in the 
case of a symmetric matrix) can never grow 
beyond that of the first vector bo . Hence we 
need not program for overflow. If the 
capacity of the machine is, let us say, 10 
decimal digits, we can fill up the cells con- 
nected with bo with 9 decimal random num- 
bers. Since the magnification factor in any 
component cannot go beyond the factor 
Vn, the provision of one empty digit in 
front of the number suffices to take care of 
the overflow, no matter how far we extend 
the procedure. 

5. Generation of the vectors b,. If the 
original matrix A is an arbitrary matrix, 
except for the fact that the reality of the 
eigenvalues is known in advance, we can 
estimate the absolutely largest eigenvalue by 
Gersgorin’s method. Taking the sum of the 
absolute values of each row (or column) and 
choosing the largest of these n positive 
numbers, we obtain an upper bound, s, on 
the eigenvalues of A. If we now put 


(5.1) 


we can be sure that the eigenvalues of this 
matrix will be limited by —2, +2. 

If, moreover, we know in advance that 
all the eigenvalues of A are positive, we put 


(5.2) C =2/] - o%. 
8 


Gersgorin’s method generally overestimates 
the largest eigenvalue of A. This has the 
disadvantage that the A-spectrum will be 
squeezed into a too small portion of the 
semicircle, leaving the remaining part un- 
used. It is advisable therefore to go through 
a preliminary Fourier analysis for a closer 
evaluation of the smallest and the largest 
eigenvalue of A, but with no demands of 
precision. For this purpose a relatively small 
value of N, e.g. 40, is sufficient, since we do 
not mind if the \-spectrum is crowded. We 
are interested only in the lower and the 
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upper limits \, and Ax (in rough approxima- 
tion). Then we put 

z as “+ 4 

(5.3) C a es Ms "Ms Me 
and thus ascertain that the eigenvalues will 
spread over the entire semicircle. 

The generation of the vectors 6, occurs 
with the help of the simple iteration scheme 
(3.5). 

6. Fourier analysis of the bh, . We isolate 
one component (or possibly two compo- 
nents) of each vector and subject this one- 
dimensional sequence of numbers y, to a 
Fourier analysis. The very first and the very 
last element, i.e. the components belonging 
to bp and by , are immediately divided by 2. 
In our subsequent formulae we will assume 
that yo and yy refer to the halved ordinates. 
The Fourier analysis consists in multiplying 
the given sequence y by the coefficients of 
a preassigned matrix I’, defined by 


(6.1) Yix = cos tk x 


N 
(6.2) “ue = ty Vka Ya 
a=0 
We can generate the coefficients y. of Ya 
concurrently with the summation, on the 
basis of the following recurrence scheme: 


yo = 1 
Y1 = Pk; 
(6.3) v2 = 2piyi — Yo 
¥s = 2prv2 — 1 
Yn = 2peyn-1 — Yn-2 
where the keyvalues p, = cos ke can be 


taken from a trigonometric table. But it is 
still simpler to let the machine evaluate and 
store the p, , by taking only the single value 
pi = cos + from the trigonometric table 
and then generate the remaining », on the 
basis of the recurrence relation (6.3). 
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We have now obtained the amplitudes 


(6.4) 


Uo, U1, te, *** , UN 


which will now be examined for peaks, fol- 
lowing the procedure of sec. 2. The presence 
of a peak is demonstrated by a + +, or 
— — sequence which interrupts the regular 
+ sequences. Such pairs are underlined 
and the calculation of the eigenvalues is 
based on these pairs, augmented by one 
neighbor to the left and to the right: 


(6.5) 


Um—1, Um, Um+i » Um4+2 - 





We form the ratio 


—_ Um—1 + 2m 1 Um+s 


(6.6) Um + Lump + Um42 
and 

ae 2-—q 

6. Be 
(6.7) °" T+¥¢ 
Then 

1 — cos (m + e) 7 

eee Sbrecars” eae 


2 


The corresponding eigenvalue of the original 
matrix A finally becomes 


(6.9) r = ys + (A be ds)’ 


Since the customary trigonometric tables 
divide the semicircle into 180°, it is conve- 
nient to choose for N, the tota] number of 
iterations, some multiple of 180. If for 
example N = 180, the angle 6; , determined 
by the integer m; plus the correction ¢;, 
is directly given in degrees. If N = 1800, 
the circle becomes divided into tenths of 


degrees (i.e. 6’), and as * gives 6; di- 


rectly in degrees. 

7. The problem of noise. We would assume 
that an iterative scheme, if pushed too far, 
might lead to a dangerous accumulation of 
rounding errors. In actual fact the genera- 
tion of the vectors , seems to be remarkably 
free of rounding errors. The rounding errors 
accumulate rather slowly and in the same 
ratio as the peaks increase in size. Hence the 
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“signal to noise ratio” does not deteriorate 
with increasing NV. Many hundreds, and even 
thousands of iterations can be performed 
without affecting the precision with which 
the 6; are obtainable. The detection of cer- 
tain very small eigenvalues demands oc- 
casionally a very large number of iterations. 
Our experience seems to indicate that the 
rounding errors have no vitiating influence 
on the result, even under such unfavorable 
conditions. 

The problem of noise can best be studied 
with the help of a matrix C which is so con- 
structed that all its eigenvalues are of the 
form 


T 
A; = 2 cos m; — 


(7.1) v 


where m; is an integer between 0 and NV 
(ef. (2.8)). In this case the field of the % 
should consist of n isolated amplitudes only, 
separated by completely blank spaces. Any- 
thing found in the blank spaces is due to 
noise. Experiments with matrices of low 
order have shown that even after hundreds of 
iterations the blank spaces remained blank 
with an accuracy of 10°. Experiments with 
large matrices and thousands of iterations 
are not yet available. 

8. Eigenvectors. The effect of the Fourier 
analysis is that of isolating one particular 
periodic component out of a mixture of 
periodic components. This is precisely our 
aim when attempting to obtain the eigen- 
vectors w; of the matrix A. As the equation 
(3.7) shows, the vector }, is a superposition 
of periodic components with amplitudes 
which are proportional to w;. This makes 
it possible to extend the technique of the 
Fourier analysis from the determination of 
the eigenvalues to the determination of the 
eigenvectors. Let a peak be near to a certain 
k = m. Then the operation 


Um = 146 bo + (cos m *) b, 
(8.1) + (cos 2m :) bo + -:> 


+14 (cos Nm 4 by 
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will generate a vector which strongly empha- 
sizes the eigenvector w; . We may thus apply 
the process of (6.2), now extended over the 
entire vector and not merely over one ele- 
ment of the vector. The vectors b, must be 
preserved in their entirety for this computa- 
tion. In contrast to the previous Fourier 
analysis, which demanded a _ systematic 
scanning for all values of k between 0 and 
N, the present process has to be carried out 
only for those specific m-values which be- 
long to the eigenvalues. 

Again the purification can be enhanced if 
we operate with second differences which in 
actual fact become second sums. In addition 
to the um-value which is nearest to a peak, we 
also evaluate the u for the left and right 
neighbors k = m — 1 and k = m + 1 and 
then take the combination 


(8.2) Um = Um. + 2Um + Umai 


This vector will be proportional to w; with 
smaller contamination on the part of the 
other w;, than if u,, alone were taken. 

9. Multiple eigenvalues. If we count the 
number of maxima of the Fourier analysis 
and find that their total number adds up to 
the order of the matrix, we may conclude 
that all the eigenvalues of our problem have 
become evident. If the total number is less 
than n, then a further investigation is needed 
concerning the missing eigenvalues. One 
possible reason for the lacunae is that the 
trial vector bp possesses accidental “blind 
spots” in the direction of some of the eigen- 
vectors. This possibility is considerably 
diminished by the fact that the components 
of bp were chosen as 9-digit random numbers. 
At all events, if we have isolated a second 
set of b, components, it is advisable to repeat 
the Fourier analysis for this set and see 
whether more peaks can be recorded in the 
second analysis. 

If this is not the case, the question of the 
missing eigenvalues can be further clarified 
by repeating the entire procedure with a 
new random vector b;. If this analysis re- 
veals new peaks, this demonstrates that our 
original trial vector had accidental de- 
generacies which were removed in the new 
analysis. But if the deficiencies still exist, we 
may assume that some of the eigenvalues are 
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multiple eigenvalues or eigenvalues which 
fall too close together. 

In this case it becomes imperative to find 
out which of the eigenvalues can be trusted 
to be single and which cannot. The decision 
can be made by the following procedure. We 
assume that we possess the Fourier analysis 
of two elements of the original trial vector, 
and the same analysis of the corresponding 
elements of the new trial vector. We form 
the second sum in the neighborhood of a 
peak, for both elements in the original set 
and compute the ratio 
us + Qui? + wit 


(9.1) r= — — 
uo, + Qu? + uG, 





Now we do the same for the new trial vec- 
tor. If r turns out to be the same in both 
cases, we have demonstrated that the eigen- 
value in question is single. If the two ratios 
are not the same, we have demonstrated that 
the eigenvalue in question represents two 
(or more) collapsing eigenvalues, or at least. 
two (or possibly more) maxima which are 
very nearly the same and which have to be 
separated by further efforts. 

10. Separation of two very close eigen- 
values. The separation of several excessively 
close eigenvalues is not an easy problem 
since two such eigenvalues operate also 
physically together and exceedingly refined 
observations (“fine structure analysis’’) 
are demanded for their separation. However, 
two very close eigenvalues can still be sepa- 
rated with relatively simple tools. We start 
with uv» and Um4: knowing that a maximum 
exists between them. However, we have 
good reasons to suspect that not one but two 
peaks can be found between these two ampli- 
tudes. We will now operate with five succes- 
sive amplitudes, viz. 

COUaD “Mao, at: We Mkts Maen 

and assume that two eigenvalues exist in 
this region. Hence we now have two p- 
values, namely p; = m + g andp = m+ e. 
We introduce 


€1€2 


II 
> 


(10.2) 


ate=eo 
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and obtain two linear equations for the de- 
termination of p and o. These two equations 
are as follows: 


(Um—2 + BUm—1 + 3m + Um4s) 
+ (QuUm2 + 3tUm—1 — Um41) 
+ (4tm2 + 3m + Uni) = 
(Um42 + 3Umsi + 3m + Um-r) 


a 
— 


(10.3) 


— (2Umi2 + 3Um41 — Um—1) 
+ (4timss + 8tUmia + Umi) = 0. 


We solve these two equations for p and o 
and then find the rdots of the quadratic 
equation 
(10.4) x —or+p=0 
which gives us «, and e. 

If the determinant of this linear set be- 
comes too small, this is an indication that 
the two peaks are too close together to be 
numerically separable (except by an increase 
of NV). In that case the two equations collapse 
into one but this one equation suffices for 
the localization of the double-peak since 
now we can put p = 0,¢ = «. 

This separation technique could be ex- 
tended to more than just two nearly equal 
peaks. The difficulty is, however, that the 
calculations become very sensitive to the 
contaminating influence on the part of the 
external peaks. For this reason we should 
first obtain the critical eigenvector by the 
Vm-method described in sec. 8. This vector is 
now a linear superposition of the two or 
more eigenvectors which belong to the 
closely bunched eigenvalues, without being 
contaminated by the external eigenvectors. 
If we use this vector as the trial vector bo 
of the recurrence scheme (3.5), the subse- 
quent Fourier analysis has a better chance 
of separating the closely spaced eigenvalues 
since the interference from the external 
peaks is removed. 

11. Very small eigenvalues. In the case 
of positive definite matrices we are some- 
times interested in the determination of a 
few excessively small eigenvalues. For very 
small eigenvalues of A (largest eigenvalue 
not exceeding 1) the relation (3.8) becomes 
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(11.1) 





_ (pry 
ms (2z) 


The appearance of N’* in the denominator 
makes even very small eigenvalues deiect- 
able. However, the case of very small eizen- 
values requires special attention, because of 
the following circumstance. The amplitudes 
u; belong to a circle and do not terminate at 
either end. The full period includes 2N terms. 
Every peak on the positive side has an ae- 
companying peak on the negative side, due 
to the relation 


(11.2) UU, = U& 
The very small eigenvalues are thus dis- 
turbed by the presence of a near peak on the 
negative side. Beyond k = 5 the influence 
from the negative side becomes negligible 
but this is not the case in the realm of the 
smallest eigenvalues. In this realm a special 
numerical scheme comes into operation. 

It is of advantage to avoid eigenvalues in 
the range between k = 0 and k = 2 alto- 
gether since in this range it is difficult to 
keep down the interference from the larger 
eigenvalues. We can push out the spectrum 
of the matrix A by putting a properly chosen 
constant in the main diagonal. We assume 
that the eigenvalues of A are already nor- 
malized to the range [0, 1]. We now modify 


A to 
one x \* 
A =A+ (x) I 


This means that the matrix C' will be defined 
by 


ye | r\ —"? 
(11.4) ¢ [2 4(%) | 4A 


The effect of this modification is that the 
eigenvalues of the new problem become 


ae «\ 
. -»+(5) 


This means—since A has no negative eigen- 
values—that the eigenvalues of the modified 
matrix cannot start below p = 2. After ob- 
taining the eigenvalues of the new matrix, 
we finally return to the original matrix by 


(11.3) 


(11.5) 
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subtracting from all \ the same constant 


(¥) 

We will now assume that the range be- 
tween p = 2 and p = 4 may contain two 
eigenvalues. The separation of these close 
eigenvalues occurs again in analogy to the 
procedure of sec. 10. Once more we introduce 
the sum and the product of the e-values, but 
with the following modification. They now 
appear in squared form: 

dd =>» 


+ € 
Moreover, «, and ¢ are not restricted to the 
range [0, 1] because the roots are measured 
from zero and we obtain for the associated 
eigenvalues directly 


(a(n) 
“i = #(5) 

Pe Tv ’ 
sl | () 


The two linear equations which determine 
co and p have coefficients which are linear 
conbinations of the first five w,-values: 
Uy, Ui, --* , Uy. The coefficients of the first 
equation are displayed in the following 
scheme: 


(11.6) 


ro 
tore 


€ 


(11.7) 


ne ne 





p —o 1 
Uo 10 
Uy 15 15 15 
(11.8) Uz 6 2 9% =0 
U3 l 9 81 


(This means that the factor of p is for 
example 


10up + 15, + bu + uz 


and so on.) The coefficients of the second 
equation become similarly: 


p —o 1 
Uo —105 
Uy 0 
Us 252 1008 4032 = 0 
Us 192 1728 15552 
Uy 45 720 11520 


If the two linear equations are not solvable 
due to the smallness of the determinant, 
we obtain only one eigenvalue but again we 
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can put p = 0, ¢ = «,, and the second equa- 
tion can be discarded. Beyond k = 4 the 
determination of « by the ratio of two second 
sums (cf. sec. 2) and the separation of two 
close eigenvalues according to sec. 10 be- 
comes applicable again. 

12. Acknewledgments. The “spectroscopic 
eigenvalue analysis” was developed in the 
winter 1953-54, during the author’s stay 
with North American Aviation, Los Angeles, 
Calif. The splendid support he received in 
his lecturing and research activities during 
this time will remain indelible in his memory. 
The almost daily discussion with Charles 
Davis and his staff was a constant source of 
inspiration. The access to the ‘‘701,” the 
electronic calculator of the I.B.M. Com- 
pany, made it possible to test the practical 
feasibility of the methods. The programming 
and coding for the 701 were planned and 
carried out by Owen Mock with great in- 
genuity. Although only matrices of low 
order (6 to 8) were employed, the generation 
of the b, vectors and the subsequent Fourier 
analysis could be studied in detail. Runs 
first of 180, then of 720, and finally of 1280, 
iterations were generated, and the non- 
accumulation of noise was demonstrated. 
Results were checked by applying another 
precision method for finding the eigenvalues 
of a matrix. 

Experiences with large-scale matrices and 
a statistical investigation of the noise prob- 
lem and its relation to the position of the 
Fourier peaks are not yet available. 
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BOTANY.—A 2-4-2 chimera of McIntosh apple. 


VoL. 45, No. 10 


Harta DerMeEn, Horticultural 


Crops Research Branch, U. 8. Department of Agriculture. 


(Received August 3, 1955) 


In 1948, a large-fruited sport of var. 
‘McIntosh’ apple, now known as ‘Kimball 
Giant McIntosh’, was reported as being a 
2-4-4 chimera (Dermen and Darrow, 1948). 
This particular chimeral makeup was deter- 
mined by examination of a shoot tip"f one 
of six one-year-old plants propagated from 
the sport portion of the mother tree. Two of 
the five remaining trees were potted and dis- 
budded to induce endogenous (adventitious') 
shoot development (for the method to induce 
endogenous shoots in the apple see Dermen, 
1948). Along with the tree that was examined 
cytologically, they were left in the green- 
house until fall when they were planted in 
the orchard. The other three trees were 
planted in the orchard when they were re- 
ceived. One of these was left intact but two 
were disbudded and died subsequently. The 
present report is based principally on the 
cytohistological study of the tree planted 
intact in the orchard. Some portions of this 
report, including illustrations of material in 
Fig. 1 have appeared previously (Dermen, 
195la, b). They are given here to elucidate 
certain points in this report and to help to 
enable the reader to follow some pertinent 
discussions. 

On one of the two greenhouse disbudded 
trees, a single endogenous bud developed and 
a cluster of shoots grew from it. On cyto- 
histological examination, the shoots were 
found to be homogeneously diploid instead 
of tetraploid as expected (Dermen, 195la). 
On the other disbudded tree, four endog- 
enous buds (a, b, c, d) developed (Fig. 1-A). 
Buds a, b, and ¢ were along a vertical line on 
the stem, whereas bud d was located at a 
point off the line from the other buds. Sub- 
sequently shoots developed from buds a and 
d only (Fig. 1-B). 

A cluster of shoots from bud a was found 
to be diploid like the cluster of shoots from 
the first disbudded tree. The shoot from bud 

1 For the use of the term endogenous in place of 


or along with adventitious, see the article by the 
author (Dermen, 1955). 


d, however, was homogeneously tetrapioid. 
These findings suggested that the young 
trees from which the endogenous shoots were 
obtained were chimerically either 2-4-2 ora 
modification of it (Dermen, 195la, b). They 
could not have been 2-44 chimeras since 
only homogeneously tetraploid shoots can 
develop endogenously from this type (Der- 
men, 1955). 

Cytohistological examination of a number 
of shoot tips from the tree grown intact in 
the orchard revealed the true condition. 
Some of these were found to have the 244 
constitution, like that previously reported 
(Dermen and Darrow, 1948); but others 
were found to be of a 2-4-2 chimeral type. A 
shoot tip of one twig with the latter consti- 
tution is illustrated in Fig. 2-A, a longisection 
through the center of the growing point. 
Layers of cells in the growing point are 
marked L-I (first layer), L-II (second layer) 
and L-III (third layer). In this material, 
L-II appears as a wide band because of its 
being composed of larger cells as compared 
with the adjacent narrow layers with smaller 
cells. Nuclei in cells of L-II were definitely 
larger than those in cells of other layers. 
This was determined by the study of see- 
tions under high magnification. It has been 
demonstrated on several occasions (Dermen, 
1947, 1951b, 1953) that large nuclei in such 
meristematic tissues are associated with 
tetraploidy. Nevertheless adequate cyto- 
logical examinations were made and the 
exact ploid nature of each apical layer and 
tissues derived from them was determined. 

In Fig. 2-B, a transverse section of the 
stem about four millimeters back of the tip 
shown in Fig. 2-A, a solid ink line was drawn 
to indicate the boundary line between tetra- 
ploid and diploid regions in the cortical tis- 
sue (for method of determination of the 
boundary line see Dermen, 1947, 1951b, 
1953). The tetraploid part of the cortex in 
this material extended from the epidermis 
to the solid ink line. The epidermis and the 
tissue inside the solid line were diploid. The 
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+ iB ie By ain a er 
Be 2 M cahi F Be it a. oe iS 
Fic. 1.—A, A disbudded one-year-old tree propagated from the large-fruited sport of a McIntosh 
apple tree, with four endogenous buds, a, b, c, and d; natural size. B, The same tree as in Fig. 


1-A photographed a month later. A cluster of shoots developed from bud a, none from 6 and c, and one 
from bud d: shoots from a were diploid and the one shoot from d was tetraploid. x. 


, 
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Fic. 2.—A, Longisection of a shoot tip from a tree also from the McIntosh sport; L-II was tetra- 


ploid and L-I and L-III were diploid. X400. 
2-A. Details are given in the text. 130. 


B, Transection of stem back of the tip shown in Fig. 
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significance of the dotted line in this figure 
is pointed out in the discussion. 


DISCUSSION 


The chimeral nature of the Kimball Giant 
McIntosh sport was first reported to be the 
24-4 type (Dermen and Darrow, 1948). 
However, the development of both diploid 
and tetraploid shoots from a disbudded tree 
propagated from the sport tree suggested 
that its constitution must have been 2-4-2 
and that the 2-4-4 chimera represented a 
form derived from the 2-4-2 chimeral con- 
dition (Dermen, 195la and b). 

It has been pointed out that endogenous 
buds on the stems of apple originated in the 
outer phloem (Dermen, 1948, 195la, 1955). 
In the stem section of the 24-2 chimeral 
twig (Fig. 2-B), the phloem, derived from 
diploid L-III, is diploid. From such a twig, 
only diploid endogenous buds can arise. If a 
shoot is constituted as a 2-4-4 chimera, only 
tetraploid endogenous buds can arise from it 
(Dermen, 1955). When both diploid and 
tetraploid shoots develop from two separate 
buds on a stem situated at different points 
of the stem circumference as shown in Fig. 1 
(buds a and d), it would indicate that in one 
part of such a stem the phloem is tetraploid 
and in another part it is diploid. This is illus- 
trated diagrammatically in Fig. 2-B. The 
dotted line in this figure was drawn to indi- 
cate that often part of the stele in the stem 
may arise from L-IJ. This was shown to 
occur in cranberry (fig. 16, 18, Dermen, 
1947) and in peach (fig. 11, 12, Dermen, 
1953). Similar developments were observed 
in apple stems but because of lack of good 
preparations for illustration, they are not 
shown. The development of diploid and 
tetraploid shoots from endogenous buds on 
the same plant is evidence that a condition 
similar to that in the cranberry and peach 
was present in the tree shown in Fig. 1. Such 
a cytochimeral makeup may have been 
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similar to that indicated in Fig. 2-B in which 
case it is assumed that growth from L-II in 
the stem had extended in part to the solid 
ink line and in part to the dotted ink line. 
Thus a shoot developed from an endogenous 
bud originating in the phloem in the area 
bounded by the dotted line would have been 
tetraploid and one from a bud in the phloem 
in the rest of the stele would have been 
diploid. 


SUMMARY 


A large-fruited sport of var. ‘McIntosh’ 
apple, ‘Kimball Giant McIntosh’ was re- 
ported earlier to have the 244 constitu- 
tion. This conclusion was based on a study 
of one tree propagated from the mother 
sport tree. When some sister trees propa- 
gated from the sport tree were used to 
produce endogenous shoots, some of the 
endogenous shoots were diploid and one 
shoot was tetraploid, indicating that origi- 
nally the sporting in the mother tree must 
have been to the 2-4-2 form. The present 
study shows this to be the case and that the 
2-4-4 condition must have represented a de- 
rived form from the 2-4-2 type. 
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MAMMALOGY.—A new subspecies of wood rat from Nayarit, Mexico, with new 


name-combinations for the Neotoma mexicana group. 


E. RaymMonp HA tt, Uni- 


versity of Kansas Museum of Natural History. 


(Received July 29, 1955) 


At lower elevations on the west coast of 
the Republic of Mexico the big, rich brown, 
dark-tailed rodent Neotoma (Hodomys) alleni 
is the common wood rat. Nevertheless, 
Neotoma mexicana also occurs there in lesser 
numbers. At the lower elevations, indi- 
viduals of N. mexicana are small. Smallness 
may be a response to warmth; anyhow 
Neotoma mexicana tends to be smaller in 
southern than in northern localities. Speci- 
mens from the coast of Nayarit pertain to 
an unnamed subspecies, which- may be 
named and described as follows: 


Neotoma mexicana eremita, n. subsp. 


Type.—Female, adult, skin with skull and 
body-skeleton, no. 64532 KU; 1 mile south of 
San Francisco, 50 feet, Nayarit; January 27, 
1955; obtained by J. R. Alcorn; original no. 
17830. 

Range.—Known only from the type locality. 

Diagnosis.—Smallest-skulled of the subspe- 
cies of Neotoma mexicana; dark grayish to dull 
ochraceous above. 

Comparisons.—From N. m. parvidens, N. m. 
eremita differs in smaller average size, less ochra- 
ceous upper parts, and less whitish (more plum- 
beous) underparts. N. m. eremita is less ochra- 
ceous even than the larger, geographically 
adajcent, N. m. tenuicauda. 

Remarks.—The holotype of this smali, dull- 
colored, rat, judged by the height of the crowns 
of the upper molariform teeth, is slightly younger 
than the holotype (71586 USNM) of parvidens 
and slightly older than the holotype (33594/ 
45629 USNM) of tenuicauda. All three are fe- 
males. Among named kinds of Neotoma, N. m. 
eremita resembles N. m. parvidens more closely 
than any other. 

Measurements.—The holotypes of eremita, 
parvidens and tenuicauda, in that order, yield 
measurements (in millimeters) as follows: 
Occipitonasal length, 39.0, 41.5, 41.7; basilar 
length, 31.6, 32.3, 33.7; zygomatic breadth, 
19.7, 20.6, ——; mastoid breadth, 15.0, 15.2, 
15.8; interorbital breadth, 4.8, 5.3, 5.4; length 


of nasals, 15.6, 15.5, 15.4; length of incisive 
foramina, 8.5, 8.6, 8.9; length of palatal bridge, 
7.0, 7.3, 7.6; alveolar length of upper molar 
series, 7.2, 7.9, 8.9; total length, 301, 295, 340; 
length of tail, 142, 141, 160; length of hind foot, 
30, 31, 31. 

Specimens examined.—Two from the type lo- 
cality. 


With assistance from the National Science 
Foundation, the Kansas University Endow- 
ment Association, and Alford J. Robinson, 
the Museum of Natural History of the 
University of Kansas has accumulated 
specimens of related kinds of wood rats, for 
example of the nominal species Neotoma 
distincta Bangs, Neotoma navus Merriam, 
and Neotoma torquata Ward as well as speci- 
mens of most subspecies of Neotoma fer- 
ruginea Tomes and WNeotoma mexicana 
Baird. Examination of these specimens and 
also of those in the United States National 
Museum including those of Neotoma trepi- 
calis Goldman gives basis for arranging all 
those mentioned above in this paragraph as 
subspecies of one species for which the oldest 
available name is Neotoma mexicana Baird 
1855. 

For example, specimen no. 63079 KU, 
here referred to parvidens, from 1 mile 
NNW of Soledad (approximately 30 km 
north of Punto Angel), 4,700 feet, Oaxaca, is 
structurally as well as geographically inter- 
mediate between N. parvidens and N. f. 
isthmica Goldman and is regarded as an 
intergrade. Of 11 specimens examined of 
N. navus from southeastern Coahuila some 
have tails as short as N. m. inornata Gold- 
man, the kind next adjacent to the north, 
and the expansion posteriorly of the frontals 
(not conspicuous in all specimens) occurs in 
some other subspecies of N. mexicana and 
leads to the conclusion that N. navus is only 
subspecifically different from Neotoma mexti- 
cana inornata. The differences between 
N. tropicalis and N. f. isthmica are of no 
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greater degree than are the differences be- 
tween some other pairs of subspecies that 
are known to intergrade—N. m. tenuicauda 
Merriam and N. m. madrensis Goldman for 
instance. It is supposed, therefore, that 
specimens from geographically appropriate 
localities will reveal that N. tropicalis does 
intergrade with M. f. isthmica. In fact, 
Dalquest (Journ. Washington Acad. Sci. 41: 
363. Nov. 14, 1951) indicated subspecific 
status under Neotoma ferruginea for Neotoma 
distincta and Neotoma torquata. Hooper 
(Occ. Pap. Mus. Zool. Univ. Michigan no. 
565: 22. Mar. 31, 1955) concluded that 
Neotoma mexicana and Neotoma ferruginea 
were conspecific because he found intergrada- 
tion between the two in specimens that he 
assigned to the subspecies fenuicauda and for 
which he used the name-combination Neo- 
toma mexicana tenuicauda. He did not, how- 
ever, implement his conclusion with name- 
combinations for the kinds (other than 
tenuicauda) affected. Doing this and doing 
the same thing for the other nominal species 
that study reveals are only subspecies re- 
duces the number of species in what has been 
referred to as the mexicana-group from eight 
to two, namely to Neotoma mexicana and 
Neotoma chrysomelas J. A. Allen. I know 
of no intergrades between these last two. 

Two other species thought by some 
students to belong in the Neotoma mexicana 
group are Neotoma goldmani Merriam, 1903, 
and Neotoma angustapalata Baker, 1951. Be- 
cause of the resemblance of their bacula to 
those of Neotoma albigula it seems best in 
the present state of knowledge to assign 
N. goldmani and N. angustapalata to the 
albigula-group. Information on the baculum 
of N.. goldmani has been published by Rainey 
and Baker (Univ. Kansas Publ., Mus. Nat. 
Hist., 7: 623. June 10, 1955). The cleaned 
bacula of N. angustapalata that Keith R. 
Kelson has shown me are indistinguishable 
from those of Neotoma albigula, figured by 
Burt and Barkalow (Journ. Mamm. 23: 
291. Aug. 13, 1942). 

The named kinds of the mexicana group 
of wood rats should stand as given below. 


Neotoma mexicana atrata Burt. 


1939. Neotoma mexicana atrata Burt, Occ. 
Pap. Mus. Zool. Univ. Michigan no. 400: 1, 
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March 1, type trom 4 miles west of Carrizozo, 
Lincoln County, N. Mex. Known only from the 
type locality. 


Neotoma mexicana bullata Merriam. 


1894. Neotoma mexicana bullata Merriam, 
Proc. Biol. Soc. Washington 9: 122, July 2, type 
from Santa Catalina Mountains, Pima County, 
Ariz. Known only from the type locality. 


Neotoma mexicana chamula Goldman. 


1909. Neotoma ferruginea chamula Goldman, 
Proc. Biol. Soc. Washington 22: 141, June 25, 
type from mountains near San Cristébal, 8,400 
feet, Chiapas, México. 

Marginal records—Chiapas: type locality. 
Guatemala: Hda. Chancol, about 13 miles 
north of Huehuetenango (Goldman, 1910: 69); 
Voleén Santa Marfa (Goodwin, 1934: 54). 


Neotoma mexicana distincta Bangs. 


1903. Neotoma distincta Bangs, Proc. Biol. 
Soc. Washington 16: 89, June 25, type from 
Texolo [Teocelo, near Jalapa], Veracruz, Mexico. 
Known only from the type locality. 

1951. N[eotoma]. flerruginea]. distincta, Dal- 
quest, Journ. Washington Acad. Sci. 41: 363, 
November 14. 


Neotoma mexicana eremita Hall. 


1955. Neotoma mexicana eremita Hall, present 
paper, type from 1 mile south of San Francisco, 
50 feet, Nayarit, México. Known only from the 
type locality. 


Neotoma mexicana fallax Merriam. 


1894. Neotoma fallax Merriam, Proc. Biol. 
Soc. Washington 9: 123, July 2, type from Gold 
Hill, Boulder County, Colo. 

1910. Neotoma mexicana fallax, Goldman, 
North Amer. Fauna no. 31: 56, October 19. 

Marginal records (Goldman, 1910: 57-58, 
unless otherwise noted).—Colorado: 5 miles 
southwest of Fort Collins; 344 miles west of 
Loveland, 5030 feet (26762 KU); Colorado 
Springs. New Mexico: Coyote Creek, 8 miles 
north of Guadalupita; Santa Rosa; Capitan 
Mountains; 1 mile south of Ruidoso, 6,500 feet. 
(35183 KU); east side of Sierra Blanca Moun- 
tains (Dice, 1942: 206); Manzano Mountains 
(V. Bailey, 1932: 183); San Mateo Peak; Datil 
Mountains (V. Bailey, 1932: 184); Grant; 18 
miles north, 1 mile east of Farmington, 6,000 
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feet (34817 KU). Colorado: Cortez (Cary, 
1911: 117); 4 miles west, 2 miles south of Cahone, 
7,000 feet (34792 KU); Bedrock (Warren, 
1942: 214); 1 mile southwest of Gateway, 4,600 
feet (34773 KU); 2% miles south of Fruita, 
4,600 feet (34772 KU); Salida (Cary, 1911: 117); 
type locality. 


Neotoma mexicana ferruginea Tomes. 


1862. Neotoma ferruginea Tomes, Proc. Zool. 
Soc. London, 1861, pt. 3, p. 282, April, type from 
Duefias, Sacatepequez, Guatemala. 

Marginal records.—Guatemala: San 
(Goodwin, 1934: 54); type locality. 


Lucas 


Neotoma mexicana griseoventer Dalquest. 


1951. Neotoma ferruginea griseoventer Dal- 
quest, Journ. Washington Acad. Sci. 41: 363, 
November 14, type from Xilitla, San Luis 
Potosi, México. 

Marginal records.—San Luis Potosf: El Salto 
(Dalquest, 1951: 363); type locality. 


Neotoma mexicana inopinata Goldman. 


1933. Neotoma mexicana inopinata Goldman, 
Journ. Washington Acad. Sci. 23: 471, October 
15, type from Chuska Mountains, 8,800 feet, 
San Juan County, N. Mex. 

Marginal records—Utah: mouth Nigger Bill 
Canyon, east side of Colorado River, 4 miles 
above Moab Bridge, 3,995 feet (Durrant, 1952: 
337). Colorado: Ashbaugh Ranch (Goldman, 
1933: 472). New Mexico: Gallup (Goldman, 
1933: 472); Zui Mountains (ibid.). Arizona: 
Canyon del Muerto and head Spruce Creek, 
9,000 feet, Tunitcha Mountains (ibid.). Utah: 
War God Spring, Navajo Mountain (Durrant, 
1952: 337). 


Neotoma mexicana inornata Goldman. 


1938. Neotoma mexicana inornata Goldman, 
Proc. Biol. Soc. Washington 61: 60, March 18, 
type from Carmen Mountains, 6,100 feet, Coa- 
huila, México. Known only from the type 
locality. 


Neotoma mexicana isthmica Goldman. 


1904. Neotoma isthmica Goldman, Proc. 
Biol. Soc. Washington 17: 80, March 21, type 
from Huilotepec, 100 feet, 8 miles south of 
Tehuantepec, Oaxaca, México. 

1910. Neotoma ferruginea isthmica Goldman, 
North Amer. Fauna no. 31: 71, October 19. 
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Marginal records (Goldman, 1910: 72).— 
Oaxaca: Coixtlahuaca. Chiapas: Teopisca; Can- 
job. Oaxaca: Juchitan; Puerto Angel; Oax:ica, 


Neotoma mexicana madrensis Goldman. 


1905. Neotoma mexicana madrensis Goldman, 
Proc. Biol. Soc. Washington 18: 31, February 2, 
type from Sierra Madre, 7,000 feet, near Gu:da- 
lupe y Calvo, Chihuahua, México. 

Marginal records (Goldman, 1910: 60, unless 
otherwise noted).—Chihuahua: Colonia Garcia; 
Pacheco; near Parral. Zacatecas: 4 miles east of 
Calabazal (Hooper, 1955: 22); Sierra de Val- 
paraiso. Durango: El Salto; 144 mile south of 
Revolcaderos (Hooper, 1955: 22); 114 miles 
west of San Luis (ibid.). Chihuahua: type locality; 
Carimechi (Burt and Hooper, 1941: 7). 


Neotoma mexicana mexicana Baird. 


1855. Neotoma mexicana Baird, Proc. Acad. 
Nat. Sci. Philadelphia 7: 333, April, type from 
{mts.] near Chihuahua, Chihuahua, México. 

Marginal records (Goldman, 1910: 56, unless 
otherwise noted).—Arizona: Rincon Mountains; 
Cave Creek, Chiricahua Mountains (Cahalane, 
1939: 435). New Mexico: Animas Mountains 
(ibid.); Guadalupe Mountains (V._ Bailey, 
1932: 182). Texas: Fort Davis; Paisano. Chi- 
huahua: Santa Eulalia. Durango: La Cienega de 
las Vacas. Sonora (Burt, 1938: 65): Cuchita; 
Oposura. 


Neotoma mexicana navus Merriam. 


1903. Neotoma navus Merriam, Proc. Biol. 
Soe. Washington 16: 47, March 19, type from 

~ Sierra Guadalupe, Coahuila, México. 

Marginal records.—Coahuila: type locality; 
Diamante Pass, 8,500 feet, Sierra Guadalupe, 
18 miles east, 3 miles north of Saltillo (36346 
KU). 

Neotoma mexicana ochracea Goldman. 


1905. Neotoma ferruginea ochracea Goldman, 
Proc. Biol. Soc. Washington 18: 30, February 2, 
type from Atemajac, 4,000 feet, near Guadala- 
jara, Jalisco, México. Known only from the 
type locality. 


Neotoma mexicana parvidens Goldman. 


1904.. Neotoma parvidens Goldman, Proce. 
Biol. Soc. Washington 17: 81, March 21, type 
from Juquila, 5,000 feet, Oaxaca, México. 

Marginal records.—Oaxaca: type locality: 
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1 mile NNW of Soledad [approx. 30 kilometers Neotoma mexicana solitaria Goldman. 


north of Punta Angel], 4,700 feet (63079 KU). 


Neotoma mexicana picta Goldman. 


1904. Neotoma picta Goldman, Proc. Biol. 
Soc. Washington 17: 79, March 21, type from 
mountains near Chilpancingo, 10,000 feet, 
Guerrero, México. 

1910. Neotoma ferruginea picta Goldman, 
North Amer. Fauna no. 31: 72, October 19. 

Marginal records (Goldman, 1910: 73).— 
Guerrero: Omilteme; type locality. Oaxaca: 
mountains near Ozolotepec. 


Neotoma mexicana pinetorum Merriam. 


1893. Neotoma pinetorum Merriam, Proc. 
Biol. Soc. Washington 8: 111, July 31, type from 
San Francisco Mountain, Coconino County, 
Ariz. 

1910. Neotoma mexicana pinetorum, Goldman, 
North Amer. Fauna no. 31: 58, October 19. 

Marginal records (Goldman, 1910: 59, unless 
otherwise noted).—Arizona: south rim, Grand 
Canyon (V. Bailey, 1935: 16); type locality; 
Winona (226370 BS). New Mexico: 10 miles 
southwest of Quemado; 10 miles west of Chlo- 
ride; Kingston. Arizona: Prieto Plateau (205483 
BS); Bradshaw Mountains (215809 BS); Sim- 
mons (215529 BS); Pine Spring (202570 BS). 


Neotoma mexicana scopulorum Finley. 


1953. Neotoma mexicana scopulorum Finley, 
Univ. Kansas Publ., Mus. Nat. Hist. 6: 529, 
August 15, type from 3 miles northwest of Hig- 
bee, 4,300 feet, Otero County, Colo. 

Marginal records (Finley, 1953: 534).—Colo- 
rado: type locality; Two Buttes [the peak]. 
Oklahoma: Tesequite Canyon. New Mexico: 
Clayton; Sierra Grande; Oak Cafion, Raton 
Range. Colorado: Trinidad; 20 miles east of 
Walsenburg. 


A. Allen. 


1898. Neotoma sinaloae J. A. Allen, Bull. 
Amer. Mus. Nat. Hist. 10: 149, April 12, type 
from Tatameles, Sinaloa, México. 

1910. Neotoma mexicana sinaloae, Goldman, 
North Amer. Fauna no. 31: 60, October 19. 

Marginal records.—Sonora (Burt, 1938: 65): 
San Javier; Mira Sol. Durango: Chacala (Gold- 
man, 1910: 61). Sinaloa: Chele (Hooper, 1955: 
22); Mazatlan (Goldman, 1910: 61). Sonora: 
Camoa (Burt, 1938: 65). 


Neotoma mexicana sinaloae J. 


1905. Neotoma ferruginea solitaria Goldman, 
Proc. Biol. Soc. Washington 18: 31, February 2, 
type from Nentén, 3,500 feet, Guatemala. 

Marginal records.—Guatemala: type locality; 
Sacapulas (Goodwin, 1942: 170-171). Honduras: 
Cerro Puca (ibid.). 


Neotoma mexicana tenuicauda Merriam. 


1892. Neotoma tenuicauda Merriam, Proc. 
Biol. Soc. Washington 7: 169, September 29, 
type from north slope El Nevada de Colima, 
12,000 feet, Jalisco, México. 

1955. Neotoma mexicana tenuicauda, Hooper, 
Occ. Pap. Mus. Zool. Univ. Michigan no. 565: 
22, March 31. 

Marginal records (Goldman, 1910: 74, unless 
otherwise noted).—Jalisco: mountains about 10 
miles north of Bolafios. Zacatecas: Plateado. 
Aguascalientes: Sierra Fria (Hooper, 1955: 22). 
Jalisco: 2 miles NNW of Magdalena (ibid.). 
Michoacan: Zamora; 9 miles southeast of Pdtz- 
cuaro, 8,000 feet (Hall and Villa, 1949: 467): 
Tancftaro, 7,850 feet (ibid.). Jalisco: type local- 
ity; Talpa; San Sebastian. 


Neotoma mexicana torquata Ward. 


1891. Neotoma torquata Ward, Amer. Nat. 
25: 160, February, type from abandoned mine 


between Tetela del Volcéin and Zacualpan, 
Morelos, México. 
1894. Neotoma fulviventer Merriam, Proc. 


Soc. Washington 9: 121, July 2, type from Toluca 
Valley, México, México. 

i894. Neotoma orizabae Merriam, Proc. Biol. 
Soc. Washington 9: 122, July 2, type from Vol- 
cin de Orizaba, Puebla, México. 

1951. Nleotoma]. flerruginea]. torquata, Dal- 
quest, Journ. Washington Acad. Sci. 41: 363, 
November 14. 

Marginel records. — Hidalgo: Encarnacién 
(Goldman, 1910: 64). Veracruz: 3 km east of 
Las Vigas, 8,000 feet (19840 KU). Puebla: 
Tehuacén, 1,700 m (Hooper, 1947: 55). 
Morelos: 2 km south of Jonacatepec (Davis 
and Russell, 1954: 77). México: north slope of 
Voleén Toluca (Goldman, 1910: 64). 


Neotoma mexicana tropicalis Goldman. 


1904. Neotoma tropicalis Goldman, Proc. 
Biol. Soc. Washington 17: 81, March 21, type 
from Totontepec, 6,500 feet, Oaxaca. Known 
only from the type locality. 
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Neotoma mexicana vulcani Sanborn. Neotoma chrysomelas J. A. Allen. 


1935. Neotoma ferruginea vulcani Sanborn, 1908. Neotoma chrysomelas J. A. Allen, Bull. 
Field Mus. Nat. Hist. Publ. 340, zool. ser., Amer. Mus. Nat. Hist. 24: 653, October 13, type 
20 (11): 84, May 15, type from south slope of from Matagalpa, Matagalpa, Nicaragua. 
Voleén Tajumulco, San Marcos, 13,200 feet, Marginal records.—Honduras (Goodwin, 1942; 
Guatemala. Known only from the type locality. 171): Montana Vasquez; Hatillo. Nicaragua; 

type locality. 
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Fia. 1.—Geographie distribution of Neotoma mexicana and Neotoma chrysomelas 


GUIDE TO KINDS: 9. Neotoma m. inopinata Goldman. . Neotoma m. pinetorum Merriam. 
1. Neotoma m. atrata Burt. 10. Neotoma m. inornata Goldman. . Neotoma m. scopulorum Finley. 
. Neotoma m. bullata Merriam. 11. Neotoma m. isthmica Goldman. . Neotoma m. sinaloae J. A. Allen. 
3. Neotoma m. chamula Goldman. 12. Neotoma m. madrensis Goldman. . Neotoma m. solitaria Goldman. 
Neotoma m. distincta Bangs. 13. Neotoma m. mexicana Baird. . Neotoma m. tenuicauda Merriam. 
. Neotoma m. eremita Hall. 14. Neotoma m. navus Merriam. . Neotoma m. torquata Ward. *& 
j. Neotoma m. fallax Merriam. 15. Neotoma m. ochracea Goldman. . Neotoma m. tropicalis Goldman. 
. Neotoma m. ferruginea Tomes. 16. Neotoma m. parvidens Goldman. . Neotoma m. vulcani Sanborn. 
. Neotoma m. griseoventer Dalquest. 17. Neotoma m. picta Goldman . Neotoma chrysomelas J. A. Allen. 











